Excess free energy and Casimir forces in systems with long-range interactions of 
van-der-Waals type: General considerations and exact spherical-model results 



Daniel Dantchev,^- * H. W. Diehl,^ and Daniel Griineberg^' 'f 

^ Fachbereich Physik, Universitdt Duisburg-Essen, Campus Essen, 45117 Essen, Germany 
^Institute of Mechanics — Bulgarian Academy of Sciences, 
Academic Georgy Bonchev Street building 4, 1113 Sofia, Bulgaria 
(Dated: February 2, 2008) 

We consider systems confined to a d-dimensional slab of macroscopic lateral extension and finite 
thickness L that undergo a continuous bulk phase transition in the limit L — > oo and are describable 
by an 0{n) symmetrical Hamiltonian. Periodic boundary conditions are applied across the slab. 
We study the effects of long-range pair interactions whose potential decays as bx~'''^'^°'^ as a; ^ oo, 
with 2 < o" < 4 and 2 < d -I- cr < 6, on the Casimir effect at and near the bulk critical temperature 
'7c,oo, for 2 < d < 4. These interactions decay sufficiently fast to leave bulk critical exponents and 
other universal bulk quantities unchanged — i.e., they are irrelevant in the renormalization group 
(RG) sense. Yet they entail important modifications of the standard scaling behavior of the excess 
free energy and the Casimir force J-c ■ We generalize the phenomenological scaling ansatze for these 
quantities by incorporating these long-range interactions. For the scaled reduced Casimir force 
per unit cross-sectional area, we obtain the form L"^ Tc /ksT « Ho(l//^oo) + fft^ (L/^oo) + 

(Jo- Z/~'^" Ho- (Iz/Coo) . Ifere Hq, H^^, and Ho are universal scaling functions; g^^ and g^ are scaling 
fields associated with the leading corrections to scaling and those of the long-range interaction, 
respectively; u) and = a + ri ~ 2 are the associated correction-to-scaling exponents, where rj 
denotes the standard bulk correlation exponent of the system without long-range interactions; ^oo 
is the (second- moment) bulk correlation length (which itself involves corrections to scaling). The 
contribution oc g^ decays for T 7^ Tc,oo algebraically in L rather than exponentially, and hence 
becomes dominant in an appropriate regime of temperatures and L. We derive exact results for 
spherical and Gaussian models which confirm these findings. In the case d + a — 6, which includes 
that of nonretarded van-der-Waals interactions in d = 3 dimensions, the power laws of the corrections 
to scaling oc 6 of the spherical model are found to get modified by logarithms. Using general RG 
ideas, we show that these logarithmic singularities originate from the degeneracy uj = uu^ — 4 — d 
that occurs for the spherical model when d-|- cr = 6, in conjunction with the b dependence of gui- 
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I. INTRODUCTION 

When macroscopic bodies are immersed into a 
medium, the forces acting between them in its absence 
are usually altered. Moreover, additional (effective) 
forces not present without the medium may be induced 
by fluctuations occurring in it. A well-known exam- 
ple of fluctuation-induced forces is the so-called Casimir 
force between metallic bodies, named after its discoverer 
H. B. G. Casimir [1], that is induced by vacuum fluctua- 
tions of the electromagnetic field and was recently verified 
through high-precision experiments [2, 3]. 

Although the Casimir effect was well received at the 
time of its discovery, interest in it diminished soon af- 
terwards, and for a long time it did not attract much 
attention. Since approximately 1970 there has been a 
resurge of interest in it, which has evolved into an enor- 
mous research activity during the past decades [4-13]. 
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There are a number of good reasons for this devel- 
opment. To begin with, fluctuation-induced forces are 
ubiquitous in nature. Casimir's original work [1] was 
concerned with the force induced by vacuum fluctua- 
tions of the electromagnetic field. Subsequently it has 
been realized that analogous forces exist that are not 
mediated by massless particles such as photons, but are 
induced by low-energy excitations such as spin-waves — 
or, more generally, Goldstone modes in systems with a 
spontaneously broken continuous symmetry — or thermal 
fluctuations. Since Goldstone modes are massless, the 
associated fluctuations are scale invariant and induce a 
long-ranged Casimir force. The same applies to thermal 
fluctuations at critical points because of the divergence 
of the correlation length. The upshot is that Casimir 
forces have turned out to be of interest for many diverse 
fields of physics, such as quantum field theories [4-8] , con- 
densed matter physics, the physics of fluids and quantum 
fluids [9-11], wetting phenomena [14-16], microfluidics, 
and nanostructured materials [17]. 

Second, owing to the progress in experimental tech- 
niques made in recent years, detailed investigations of 
Casimir forces have become possible [2, 3, 18-23]. Third, 
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a further important reason for the ongoing interest in 
Casimir forces is that they exhibit universal features: 
Microscopic details of both the fluctuating medium as 
well as the immersed macroscopic bodies do not normally 
matter, at least as long as long-range interactions are ab- 
sent or may be safely ignored. 

Last but not least, an equally important reason has 
been the theoretical progress in dealing with interact- 
ing field theories with boundaries that has been achieved 
since the 1980s [24-30]. This has led to detailed investi- 
gations of the Casimir effect for interacting field theories 
[13, 28, 31- 35]. 

In this paper we will be concerned with the thermo- 
dynamic Casimir effect — i.e., the Casimir effect induced 
by thermal fluctuations. Our aim is to study the effects 
of long-range interactions of van-der-Waals type on the 
Casimir force in systems undergoing a continuous bulk 
phase transition. To this end we shall consider long-range 
two-body interactions with a pair potential v^'^^ (x) that 
behaves as 

v'-'^^x) « const a;- ('^+'^) (1.1) 

in the large-distance limit. The familiar dispersion forces 

in fluids belong to this category: Important examples are 
the nonretardcd and retarded van-der-Waals interactions 
of a d = 3 dimensional fluid, which correspond to the 
cases cr = 3 and a = 4, respectively. 

According to scaling considerations (to be recalled in 
Sec. II), the leading infrared singularities at the criti- 
cal point of systems with short-range forces do not get 
modified by such long-range interactions because the as- 
sociated pair potentials decay sufficiently fast at large 
distances. They are irrelevant in the renormalization 
group (RG) sense, giving corrections to the leading crit- 
ical behavior. Long-range interactions of this kind have 
been termed "subleading long-range interactions" [36- 
39]. They are generically present in fluids [40-42] but 
occur also in other, for example, magnetic systems. 

For typical three-dimensional systems of the n-vector 
type with n < oo, the associated correction-to-scaling 
exponent Wa is larger than the familiar exponent u> 
that governs the leading corrections-to-scaling (see, e.g., 
Ref. [43]). Hence such long-range interactions yield next- 
to-leading corrections-to-scaling. 

Despite their irrelevance, they have important conse- 
quences, even for the near-critical behavior of bulk sys- 
tems. Since they involve pair potentials that decrease as 
inverse powers of the distance x in the limit x — > oo, the 
usual exponential large-x decay of correlations away from 
the critical point gets replaced by an algebraic one. 

Their consequences for the medium-induced force be- 
tween two macroscopic bodies immersed into the medium 
a distance L apart is of a similar kind and importance: 
They yield contributions that decay quite generally as 
an inverse power of L, irrespective of whether or not the 
temperature T is close to the bulk critical temperature 
Tc. CO of the medium. When T ~ Tc.ca, they compete 
with the long-ranged Casimir force produced by critical 



or near-critical fluctuations. As previous work [36-39, 44] 
suggests, and will be shown in detail below, they actually 
become the dominant part of the medium-induced force 
in a certain regime of temperatures and L. 

We will consider the case of a slab geometry of cross- 
sectional area A = L'^\~^ and thickness L. Reliable results 
for this geometry are important for the interpretation 
of Monte Carlo simulations of appropriate models with 
subleading long-range interactions. 

In view of our above remarks, a most obvious system 
class to consider would be fluids. To describe the long- 
distance physics of classical fluids near their liquid-gas 
critical point, a one-component order parameter is used. 
Instead of considering this case, we will focus our at- 
tention on systems that involve an n-componcnt order 
parameter and can be modeled by an 0{n) symmetrical 
Hamiltonian, and investigate them in the limit n — > oo. 
For simplicity, we will restrict ourselves to the case of pe- 
riodic boundary conditions along all — namely, both the 
perpendicular as well as the d — 1 principal parallel — 
directions. Under these conditions, the large-n limit 
of the 0{n) model is equivalent to the spherical model 
[45-47]. We will present exact results for the Casimir 
force at and above the bulk critical temperature T^, both 
for spherical and Gaussian models with subleading long- 
range interactions. 

Our motivation for considering spherical models is 
twofold. First, studying the effects of such long-range 
interactions on the Casimir force for such models is an 
interesting problem in its own right. Second, the exact re- 
sults obtained for these models provide nontrivial checks 
for the results of perturbative field-theoretic renormaliza- 
tion group approaches and are expected to give valuable 
guidance for acceptable approximations, an issue we plan 
to take up in a subsequent paper [48]. 

A special feature of the spherical model with 2 < d < 4 
is that the correction-to-scaling exponents uja- and to be- 
come equal when d -|- cr = 6, a condition satisfied, for 
example, for nonretarded van-der-Waals interactions in 
d = 3 dimensions. As our exact results show, the 
corrections-to-scaling induced by the long-range interac- 
tion (1.1) then get modified by logarithms. 

The remainder of this paper is organized as follows. In 
the next section, we provide the required background on 
Casimir forces. We begin by recalling the definition of 
the Casimir force. Then wc discuss its scaling form when 
all interactions are short ranged, specify the form of the 
subleading long-range interactions to be considered, and 
recapitulate the scaling arguments which show that they 
do not modify the leading critical singularities. Next, we 
generalize the finite-size scaling ansatz by incorporating 
them. In Sec. Ill wc introduce the spherical model with 
subleading long-range interactions which we solve for 2 < 
d < 4 to produce exact largc-n results for the Casimir 
force. The finite-size behavior of the equation of state is 
analyzed in Sec. IV. Section V deals with the finite-size 
behavior of the free energy and the Casimir force. Sec. VI 
contains a brief summary and discussion. Finally, there 
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are three appendixes in which various technical details 
are explained. 

II. BACKGROUND 
A. Definition and scaling form of Casimir force 

We consider a statistical mechanical system, a model 
magnet or fluid, whose shape is a d-dimensional slab of 
thickness L and hyperquadratic cross-section with area 
A = As previously mentioned, we choose peri- 

odic boundary conditions along all d principal hypercu- 
bic axes, so that the system has the topology of a d-torus. 
Unless stated otherwise, the dimensionality d is presumed 
to satisfy 2 < c? < 4. 

Let Fl,a{T) be the total free energy of the system. 
Taking the thermodynamic limit ^ oo at fixed L < 
00, we denote the reduced free energy per cross-sectional 
area A as = liraA^oo FL,A/AkBT. For L — > oo, 

f]^{T)/L approaches fhk{T), the reduced bulk free energy 
density [49]. We therefore introduce the reduced excess 
free energy by 

/ex(T,L)-/L(T)-L/bk(r) . (2.1) 

The limit L ^ oo of this quantity exists, but depends 
on the boundary conditions: for periodic boundary con- 
ditions and the film geometry with boundary planes 5Bi 
and *B2 introduced above, we have [49] 

U{T,H,oo) = l^/ periodic be 

^ ' ' ' \ fs,i + /s,2 , film geometry, ^ 

where fg^i, i = 1,2, are the surface excess free energy of 
the respective semi-infinite systems bounded by 5Bj. 

In either case, the thermodynamic Casimir force per 
unit area is defined in terms of /ox as 

J^c{T,L) = -kBT^I^^^^. (2.3) 

According to this definition, this quantity is a gen- 
eralized force conjugate to the thickness L of the slab, 
which approaches zero as L — > oo. We are interested 
in its behavior for L ^ a, where a is a typical micro- 
scopic length scale (we henceforth set to unity) . Suppose 
for the moment that all interactions are short-ranged. 
Then finite-size scaling theory should be applicable in 
this limit. According to it, the Casimir force takes the 
scaling form [11, 50] 

J'ciT, L)/kBT = L-" So(L/eoo) , (2.4) 

where ^oo is the bulk correlation length [51], while Sq is 
a universal scaling function. This holds up to eventual 
contributions from regular background terms and irrel- 
evant scaling fields, which we disregard for the moment 
but will come back to later, in particular, in Sec. II B 2. 



As the temperature T approaches its bulk critical value 
^c.ooi with L fixed at a finite value, the correlation length 
^oo diverges and L/^oo ^ 0. The corresponding limiting 
value of the scaling function Sq (which exists) is conven- 
tionally written as 

So(0) = (d-l)Ac, (2.5) 

which defines the so-called Casimir amplitude Ac [14]. 
This quantity is related to the critical Casimir force via 

^c(Tc,oc,i)ABTe,^ = (rf-l)^. (2.6) 

Just as the scaling function Sq, it is a universal quantity; 
it is independent of microscopic details, but depends on 
the bulk universality class considered and on other gross 
features such as boundary conditions. 

Let us be a bit more precise. Suppose that instead 
of choosing periodic boundary conditions we considered 
a lattice model with free boundary conditions along the 
perpendicular direction. Then the topmost and lowest 
layers of the system would be free surfaces, correspond- 
ing to macroscopic planar boundaries between which the 
Casimir force acts. Provided (a) no symmetry-breaking 
boundary terms are included in the Hamiltonian and (b) 
no long-range surface order is possible for T > Tc^ooi 
one expects the long-distance physics of the system near 
the bulk critical point to be described by an 0{n) (j)'^ 
model with Dirichlet boimdary conditions. This is be- 
cause upon coarse graining, the lattice model with free 
boundary conditions maps onto such a continuum field 
theory, albeit one satisfying Robin boundary conditions 
inside of averages [29, 30]. 

If conditions (a) and (b) are satisfied, one has reason 
to believe that the theory belongs to the basin of attrac- 
tion of the fixed point describing the so-called ordinary 
surface transition. This fixed point is infrared-stable and 
corresponds to a Dirichlet boundary condition on large 
scales. The analogs of the Casimir amplitude Ac and 
the scaling function So for this case of Dirichlet bound- 
ary conditions on both surface planes differ from their 
counterparts for periodic boundary conditions. Details 
of the mesoscopic Robin boundary condition or micro- 
scopic details of the boundaries do not matter as long 
as the resulting continuum theory belongs to the basin 
of attraction of the mentioned fixed point. 

More generally, we have for a film geometry bounded 
in one direction by a pair of parallel boundary planes *Bi 
and 5B2 the following situation. Universal quantities such 
as the Casimir amplitude Ac or the scaling function Sq 
depend (for given bulk universality class and short-range 
interactions) on gross properties of both boundary planes. 
Let SUCj denote the universality class pertaining to the 
surface critical behavior of the semi-infinite system with 
boundary plane 03 j ("surface universahty class" SUC), 
where i = 1 or 2. To specify universal quantities like 
the Casimir amplitude, we can write ^^uCi,suC2^ where 
possible choices of SUCi and SUC2 are "ord", "sp", and 
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"norm" , the SUC of the ordinary, special, and normal (or 
extraordinary [52]) transition, respectively. The above- 
mentioned case of Dirichlet boundary conditions on *Bi 
and 5B2 corresponds to the choices SUCi = SUC2 = ord. 

Systems with 0(n)-symmetrical Hamiltonian and 
short-range interactions have been studied in such film 
geometries for various choices of SUCi and SUC2 by 
means of the e-expansion about the upper critical di- 
mension d* = 4 [9, 31 35], Monte Carlo simulations 
[10, 53, 54], and other techniques [11]. A fairly up-to- 
date survey of pertinent results may be found in the latter 
reference. More recent results arc contained in Rcf. [54]. 
Aside from these cases and the one of periodic boundary 
conditions, also slabs with aniiperiodic boundary condi- 
tions have been considered for systems with short-range 
interactions [31, 32]. 

Going back to the case of periodic boundary condi- 
tions, we now turn to the question of how to include 
subleading long-range interactions. 

B. Subleading long-range interactions 

We consider long-range two-body interactions with a 
pair potential v^'^\x) of the kind (1.1). Let us begin by 
recalling how the relevance or irrelevance of such inter- 
actions for bulk critical behavior can be assessed. 



1. Relevance/Irrelevance criterion 

Let Hsr be the standard <j)'^ Hamiltonian represent- 
ing the bulk universality class of the n-vector model 

with short-range interactions for d below d* — 4, its 
upper critical dimension. At the bulk critical point, 
the n-component order parameter field transforms as 
(p £"^["^10 under changes /i fi£ of the momentum 
scale, where the scaling dimension A[(/i] is given by 

A[0] = (d-2 + ?7)/2 . (2.7) 

Adding to Tisr a long-range interaction term with pair 
potential v^^^x), we consider the Hamiltonian 

n = nsr + b j O^"^ {x) d'^x , (2.8) 

where 0^'^^{x) denotes the nonlocal operator 

0(-)(a^)= y'd'^J/^;(-)(y)0(a;-|)c/>(a; + |), (2.9) 

and h is the associated coupling constant. 

We now ask under what conditions the short-range 
fixed point remains infrared-stable with respect to this 
O^'^^ perturbation. Upon insertion of the limiting form 
(1.1) into it, we can use Eq. (2.7) to conclude that the 
scaling dimension of the associated scaling operator, at 
the short-range fixed point, is given by 

A[0('^)] = d - 2 -h r? -h c7 . (2.10) 



The corresponding scaling field ^ b varies as £ y" in 
the infrared limit ^ — > 0, with the RG eigenexponent 

y^ = -LJ^ = d- A[0^''^] = 2-7]- a . (2.11) 

Depending on whether the correction-to-scaling exponent 
Wo- > or < 0, the short-range fixed point is locally 
stable or unstable to such perturbations. Hence wc ar- 
rive at the following irrelevance/relevance criterion: The 
long-range perturbation a 6 is irrelevant at the short- 
range fixed point if 

a>2-r), (2.12) 

and relevant if a < 2 — t]. Note that here and elsewhere 
in this paper, rj always means the correlation exponent 
of the short-range case. 

The case when this criterion suggests these long-range 
interactions to be relevant has been studied in the liter- 
ature in the context of bulk critical behavior. For a < 2, 
the upper critical dimension above which Landau theory 
holds is lowered from d* = 4 to di^{(7) = 2g. In the 
regime a < d < dl^(a)^ the values of the critical expo- 
nents depend on tr, where the analog of r/ is given exactly 
by ??ir = 2 — (T [55-60] . For given d, a crossover from the 
critical behavior characterized by these critical exponents 
to one representative of systems with short-range inter- 
actions is predicted to occur at a = 2 — 77 [56, 61-64]. 
This crossover has recently been reexamined for d = 2 
by numerical means [65]. 

Since we assume in our subsequent analysis that 2 < 
(7 < 4, the irrelevance criterion (2.12) is satisfied. Associ- 
ated with the long-range interaction (2.9) therefore is an 
irrelevant scaling field go- ~ 6 whose RG eigenexponent 
is given in Eq. (2.11). We next generalize the finite-size 
scaling ansatz for the free energy by incorporating 

2. Finite-size scaling 

Allowing a magnetic field H to be present, we con- 
sider the reduced free energy per unit cross-sectional area 
A = Z/^~^ of the previously specified slab with peri- 
odic boundary conditions, in the thermodynamic limit 
L± 00. According to the phenomenological theory of 
finite-size scaling [50, 66, 67], this quantity can be decom- 
posed into a regular background contribution f^^{T, H) 
and a singular part f2^^(T,H): 

MT, H) = /f^(T, H) + f;^{T, H) . (2.13) 

This decomposition entails analogous decompositions of 
the bulk and excess free-energy densities fi,]^{T,H) and 
/ox(r,i7,i), respectively. 

Before turning to the singular parts, let us briefly com- 
ment on the regular background terms. For simple lattice 
systems with short-range interactions it has been found 
that the regular background terms of the excess free en- 
ergy in the case of periodic boundary conditions agree to 
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high accuracy with those of the bulk free energy [50]. 
This is understandable: Periodic boundary conditions 
preclude surface and edge contributions to the total free 
energy and hence terms of this kind that are analytic in 
temperature and magnetic field. Yet, it must be remem- 
bered that free energies and their regular background 
contributions are no universal properties, but depend on 
microscopic details of the system considered. Suppose a 
given system with periodic boundary conditions that be- 
longs to the bulk universality class of the d-dimensional, 
n-component (j)^ model. Then we can choose a simple 
lattice n- vector model with nearest-neighbor interactions 
to investigate its universal critical behavior. However, in- 
clusion of any irrelevant interaction — in particular, long- 
range interactions — that were dropped when making the 
transition from the original system to the lattice model 
is expected to modify the regular background contribu- 
tions of the (bulk and excess) free energy. In other words, 
the empirical fact that the regular background contribu- 
tions of the bulk and excess free energies of simple lattice 
models with short-range interactions can be chosen to 
be equal when periodic boundary conditions are applied, 
does not imply that the same is true for microscopically 
more realistic model with additional (irrelevant) inter- 
actions. In particular, this must be kept in mind when 
adding irrelevant long-range interactions. 

The singular parts /""s, /^j^"s, and should have a 
scaling form. Specifically, /|™^(T, H, L) should take the 
finite-size scaling form 

/|»g(T,/i,L) 

(2.14) 

on sufficiently large length scales, where uj is the previ- 
ously mentioned standard correction-to-scaling exponent 

of short-range systems. Further, gt, gh, gu>, and g„ denote 
scaling fields. The first two are the leading even and odd 
relevant bulk scaling fields (namely, the "thermal" and 
"magnetic" scaling fields). For simple magnetic systems 
they behave as 

gt^att, t= {T -Tc,oo)/Tc,oo, (2.15) 

and 

gh ^ ahh, h = H/kBT^oo, (2-16) 

near the bulk critical point {T,H) = (Tg^oojO), where at 
and tth are nonuniversal metric factors; for fluid systems, 
both become linear combinations of t and the devi- 
ation of the chemical potential from the critical point, 
because of "mixing" (see, e.g., Refs. [68, 69]). For sim- 
plicity, we will use magnetic language and work with the 
above expressions henceforth. 

Likewise, the previously introduced scaling field asso- 
ciated with the long-range interaction (2.9) is expected 
to vary as 

ga « a^b (2.17) 



for small b. The ellipsis in Eq. (2.14) stands for analogous 
expressions involving further scaling fields, all of which 
we assume to be irrelevant; this means, in particular, 
that all relevant scaling fields other than gt and gh are 
taken to vanish. Moreover, we assume that none of the 
suppressed irrelevant scaling fields is dangerous irrelevant 
(see, e.g., Appendix D of Ref. [70]), so that all of them 
may be safely set to zero. 

Current estimates of the correction-to-scaling expo- 
nent uj{n,d) of the d-dimensional n-vector model give 
a;(l,3) ~ 0.81 and somewhat smaller values for n = 2 
and n = 3, such as ci;(3,3) ~ 0.80 [71-73]. On the other 
hand, the well-known exact spherical-model (SM) value 
is 

(2<d<4)==w(oo,2<d<4)=4-d. (2.18) 

Let us compare these numbers with the appropriate 
analogs for the correction-to-scaling exponent one 
can derive from Eq. (2.11). The cases of nonretarded 
and retarded van-der-Waals interactions in d dimensions 
correspond to the choices a = d and a = d+1, giving 
cod = d — 2 + T] and oJd+i = d — 1 + t], respectively. 
Both exponents are positive in the regime of dimensions 
2 < (i < 4 we are concerned with. For finite n, where 
77 > 0, the latter remains larger than iv in this whole 
regime, whereas u)d would become smaller than u slightly 
below d = 3. In the spherical limit n oo, this sign 
change of — a; occurs aX d — '3 where Ud = u = 1. In 
our analysis of the spherical model given below we shall 
first assume that d + a < 6. Then the possibility that 
w > is ruled out. The borderhne case c?-|- ct = 6 of the 
spherical model is special because w = w^- = 4 — d. Ow- 
ing to this degeneracy, it requires special attention and 
will be discussed separately. 

For the time being we therefore take it for granted 
that the irrelevant scaling fields go, and g^ yield leading 
and next-to-leading corrections to scaling in the criti- 
cal regime, respectively. However, away from the bulk 
critical point, the long-range interaction is expected to 
modify the large-i behavior of /ex and the Casimir in 
a qualitative manner so that they decay as inverse pow- 
ers of L rather than exponentially [44]. To see how this 
translates into properties of the scaling function X, let 
us denote the scaling variables appearing in Eq. (2.14) as 



t 




h 










T — UJ 

= 9wL 



and expand X as 

X{i,h;ga,gu,) = Xo{t,h) + g^ X^{t,h) 

+ g^X^{i,h) + ... , (2.20) 

where it is understood that all suppressed scaling fields 
have been set to zero. 
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The scaling functions Xq and X^^ obviously arc prop- 
erties of the short-range universality class. A similar, 
though somewhat more restricted statement applies to 
Xcr: Just as the other two, it may be viewed as the ex- 
pectation value of a quantity, computed at the infrared- 
stable fixed point of the (/>^ model with short-range inter- 
actions in a periodic slab of thickness L = 1. However, 
it differs from those inasmuch as, in its case, this quan- 
tity is the nonlocal operator 0^'^^{x) associated with the 
long-range interaction, whereas the other do not involve 
this interaction at all. 

At the bulk critical point gt = gt = 0, all three of these 
scaling functions are expected to take finite, nonzero 
values. Specifically, the critical value of Xq yields the 
Casimir amplitude: 



Ac = Xo{0,0) . 
We denote its analogs for X^ and X^ as 
A^,c = ^0,(0,0), 



(2.21) 



(2.22) 
(2.23) 



The former controls the leading corrections to the asymp- 
totic behavior of the critical excess free energy, the latter 
its contribution linear in g^ originating from the long- 
range interaction (2.9). 

Next, we turn to a discussion of the behavior as L ^ oo 
when T > Tc^oo- In this limit, either the scaling vari- 
able t, or both t and /i, tend to infinity. As explained 
above, both functions Xq and must decrease as 
~ exp[-L/^(^^)(T,H)], where ^'■^'\T,H) is the true cor- 
relation length of the system with short-range interac- 
tions. Let us set gh = for the sake of simplicity. As 
t — > oo we then should have 



Xo {i, 0) _ ~ exp - I const | P O (In i) 



(2.24) 



and similar asymptotic behavior for X^,. However, for 
the function Xcr{i, 0) we anticipate the limiting form 



(2.25) 



The exponent C, introduced here characterizes the asymp- 
totic dependence on L/^ via X„ ~ [L/S)~'^. Our results 
for both the spherical {n = oo) and the Gaussian model 
(GM) derived in the following sections yield 



CsM = Cin = oo) = Cgm = 2 , 



(2.26) 



in conformity with Ref. [44]. 

In the regime L/£^^ 1 where Xq and X^ are exponen- 
tially small, the implied contribution ~ g„ to the excess 
free energy should become dominant: 

/ex(L5r » 1) « 9a L-(<i+-+'7+C-3)^-K _ (2.27) 

and imply a corresponding large-L behavior 

J^c ~ gaL-^'^+'^+^+i-^^gr^ (2.28) 



of the Casimir force. 

In the cases of the spherical and Gaussian models, 
where 77 = and ( is given by Eq. (2.26), the large- 
ly dependence of /ox reduces to ~ 2^-(<i+<^-i). Our exact 
results for the spherical and Gaussian models given below 
confirm these findings. In fact, there are reasons to ex- 
pect that the latter L dependence applies more generally 
even when t] > 0. As proven some time ago by lagolnitzer 
and Souillard [74], using the Griffiths-Sherman-Kelly in- 
equalties [75] , the two-point net correlation function of a 
ferromagnetic system whose interactions decay as v^'^^ (x) 
in Eq. (1.1) cannot decay faster than the potential. Al- 
though we are not aware of any rigorous proof that they 
cannot decay slower than the potential either [76, 77], 
it seems most natural to us to assume that this cumu- 
lant decays as x ^ 00 according to the same power law 
as the interaction potential, barring eventual logarithmic 
corrections in special cases. 

Now, the correlation function 



G{x)^{S{x)S(0))-{S{x){S{0)) 



(2.29) 



of a slab of size oo''"^ x L under periodic boundary con- 
ditions (pbc) can be expressed in terms of its bulk coun- 
terpart Goo via 



G' 



(pbc) 



i^) = Goo{x- jLei) 



(2.30) 



where ei is a unit vector along the finite 1-direction. For 
fixed X, the terms with j yield L-dependent devi- 
ations from the j = bulk term that decay ^ L~''^~'^ 
as L ^ oc', provided the large-distance behavior (1.1) of 
v^"^^ carries over to Goo- This suggests that (away from 
criticality) the excess contribution to the free energy of 
a d-dimensional volume is down by a factor L~'^, so that 
the excess density /ox behaves as L^^'^^'^. To ensure con- 
sistency with Eq. (2.27), we must therefore have 



C = 2-r,. 



(2.31) 



Although our results for the spherical and Gaussian 
models described below are in conformity with this pre- 
diction, they do not provide a nontrivial check of it be- 
cause 77 vanishes. Such a check should in principle be 
possible within the framework of the e expansion. To 
this end, one would have to compute the scaling func- 
tion X^ using RG improved perturbation theory to suf- 
ficiently high orders, verify its limiting behavior (2.25), 
and confirm its consistency with Eq. (2.31). 



III. SPHERICAL MODEL 

A. Definition of the model 

Let £ C Z'^ be the set of sites a; of a simple hyper- 
cubic lattice of size L1XL2X • • • x^d- Imposing periodic 
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boundary conditions along all d principal directions, we 
consider a spherical model with the Hamiltonian 



UbT 2 ^ UbT 



-hY,Six)+sJ2s\x) (3.1) 



whose spin variables S{x) e M satisfy the mean spherical 
constraint 



(3.2) 



Here the cardinality of the set £, is the total number 

of sites (or spins). Further, s is a real positive variable, 
called spherical field, whose value is to be determined 
from Eq. (3.2). For systems such as the one considered 
here, whose spins are all equivalent by translational in- 
variance, the constraint (3.2) fixes all averages {S^{x)), 
yx e £, to be unity. 

As before, h = H/ksT denotes a reduced magnetic 
field. The pair interaction J{x) consists of nearest- 
neighbor bonds and a long-ranged contribution of the 
type v^'^'^ specified in Eq. (1.1), with 2 < a < 4; we use 
the choice 



J{x) = Ji da:,l + 



(/9g+a;2)(<i+<^)/2 



(3.3) 



with Ji > and J2 > 0, where po > sets a crossover 
length scale beyond which J{x) varies approximately as 
h'. 



Jox-'^-". 



B. Properties of the interaction potential 

In Appendix A wo show that the Fourier transform 



J(q) = ^J(a;)e-' 



of this interaction can be written as 

J{q) = J{0)-KkBTVL{q) 

with 



K = - ^ Mg) 



(3.4) 



(3.5) 



(3.6) 



q=0 



where 0,{q) behaves as 

d 

n{q) = g2 _ bq- + hq^ + 64,1 ^ + o{q^) (3.7) 



a=l 



for small q, and > 0, 6 > 0, 64 > 0, and 64 + 64^1 > 0. 
The term oc &4_i is anisotropic in g-space. It is a conse- 
quence of the fact that the hypercubic lattice breaks the 



Euclidean symmetry down to the symmetry of a hypcr- 
cube. For other, less symmetric lattices more than two 
fourth-order invariants and hence additional anisotropic 
terms would appear. 

Owing to our choice (3.3) of interaction constants, we 
have J{x) > for all lattice displacements x. A straight- 
forward consequence is that the Hamiltonian (3.1) has a 
unique ground state whose energy for /i = is given 
by J(0). Furthermore, J(0) > J{q) for all nontrivial 
wave- vectors q in the first Brillouin BZi. It follows that 
the resulting values of 6, . . . , 64,1 must be such that the 
equation 1 — bq'^^^ + q^[b4 + 64,1 ^^^(^a/?)^] = has no 
real- valued solutions q. 

Since the nonanalytic contribution ~ q'^ arises from the 
large-distance tail of J{x), its coefficient ksTKh should 
not dependent on the details of how J{x) behaves at 
small distances and hence be independent of po- Our 
result 



ksTK 7r'^/2r(-cr/2) 



J2 



2<'T[{d + a)/2] 



(3.8) 



derived in Appendix A, confirms this expectation. On 
the other hand, the coefficients of the analytic terms of 
orders q^ and of J{q) depend, of course, on Ji and pQ. 

Note that the Fourier transform of the second term on 
the right-hand side of Eq. (3.3) yields a contribution to 
the nearest-neighbor coupling J(a;)|a;=i that depends on 
K. This dependence can be utilized to modify this con- 
tribution and hence the nearest-neighbor coupling for a 
given value of J2 by varying k. If we choose, for sim- 
plicity, the value zero for the coupling constant Ji in 
Eq. (3.3), then the parameter h becomes (cf. Appendix A) 

6 = --r(-(T/2)r(2-cr/2) sin(7r<T/2) (po/2)''"^ , (3.9) 

TT 



which reduces to 



&=gP0, fi = cr = 3, 



(3.10) 



for the case of nonretarded van-der-Waals interactions in 
three dimensions. 



C. Solution of the model, free energy, and 
constraint equation 



Defining 



^^2 = — In Jig) =— 



we introduce the parameter 



1 

V2 



2skBT 



J(0) 



and the mode sum 



(3.11) 



(3.12) 



(3.13) 



qSBZi 
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where L = (Li, . . . , Ld). As is shown in Appendix A. the 
coefficient for our choice (3.3) of interaction constants 
takes the value 



V2 



Pi 



2(cr - 2) 



(3.14) 



when Ji = 0. 

Expressed in terms of the above quantities, the total 
free energy Fl{K, h) of our model is given by [11] 

f^|>=/.».(A-) + isup{2a,„ML)-i^.-|l} 

(3.15) 

with 



, K K 

In 

27r V2 



(3.16) 



To determine the required suprcmum, wc differentiate 
Eq. (3.15) with respect to r. This yields as condition 
from which r = rL{K, h) — or, equivalently, the spherical 
field s of Eq. (3.12) — is to be determined, the constraint 
equation 



with 



K 



Krl 



- y 

10,1 ^ 



(3.17) 



(3.18) 



The latter quantity is obviously related to Ud,Q, via 
Ud,a{rL\L) = C/d,Q(0|i) + ^ ^'^^ Wd,n{x\L) dx. (3.19) 

Let us recall that the constraint equation (3.17) can 

be recast in the form of an equation of state [11]. To sec 
this, note that Eq. (3.15) yields for the magnetization 
density niL the result 



mr.{K, h) - — j^^^^j^ j^^ , 



(3.20) 



whenever the suprcmum is attained for the solution 
of Eq. (3.17). Using this to eliminate rz in favor of rriL 
and h gives us the equation of state 



(1 -ml)K = W, 



d,n 



rriLK 



In view of this correspondence between the constraint 
equation (3.17) and the equation of state (3.21), we will 
take the liberty of referring to the former henceforth as 
the equation of state. 

Prom Eq. (3.20) one can easily read off that tl for 
/i = has the familiar meaning of an inverse susceptibil- 
ity. Let us define the susceptibility by 



XL{K,h) 



dniLiK, h) 
dh 



(3.22) 



Taking the derivative of the above-mentioned equation 
with respect to ft. at /i = then gives the desired relation 



(3.23) 



We are interested in the limit where all linear dimen- 
sions L2, ■ . ■ , Ld ^ 00 while Li remains fixed at the finite 

value Li = L. Let us employ the following convenient 
convention: Whenever the bold symbol L in quantities 
such as Wd,n{r\L) or has been replaced by L, it is un- 
derstood that the so specified thermodynamic limit has 
been taken. For instance, C/d.n(r|L) stands for 



lim 



UdMr\L) 



(3.24) 



and Ud^niTLlL) means this fimction, taken at the cor- 
responding limiting value tl = liniL2,...,id^oo '^^ 
supremum r^, i.e., of the solution to Eq. (3.17). 



IV. FINITE SIZE BEHAVIOR OF THE 
EQUATION OF STATE (3.17) 

A. Decomposition of mode sums into bulk and 
size-dependent contributions 

In order to determine the finite size behavior of the 
excess free energy and its consequences for the Casimir 
force, we must investigate the L dependence of the mode 
sums Ud,Q{rL\L) and Wd,QirL\L) for large L. Let us 
first focus our attention on the explicit L dependence 
of these quantities by considering them at an arbitrary 
X-indepcndent value of r. Writing 



UdMr\L) = C/d,n(r|oo) + Af/d,n(r|L) 



and 



Wd,n{r\L) = VFd,n(r|oo) + AWd,n{r\L) , 
we split off their L-independent bulk parts 

rid) 



(4.1) 
(4.2) 



and 



t/d,n(r|oo) = i / ln[r + n{q)] (4.3) 



(3.21) where 



qeBZi 



n 



dqa_ 

27r 



(4.5) 



is a convenient short-hand, from their L-dependent rests 
AUdSA{r\L) ^■iid AWd,n{r\L). Using Poisson's summa- 
tion formula (A4) (see Appendix C), the latter can be 
written as 

AUd,n{r\L) = ^ / cos {q^kL) ln[r + n{q)] (4.6) 
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and 



respectively. 



(4.7) 



B. Bulk equation of state 

Next, we consider the equation of state (3.17) in 
the bulk limit L oo. At the bulk critical point 
K = Kc,L=<x, h = 0, its solution r = rz,=oo must van- 
ish. Hence the critical coupling K^. is given by 



(4.8) 



As indicated, this quantity depends on the interaction 
parameter b as well as on all other interaction parameters 
64, 64,1, etc appearing in fl{q). Defining the scaling fields 

9t and gh as 



9t 



Kc,oo - K , gh = h/\/K 



(4.9) 



wc find from Eq. (3.17) that the bulk quantity Too is to 

be determined from 

-9t = {Qh/roof + Wd,n(roo|oo) - Wd,a(0|oo) , (4.10) 
the "bulk equation of state" . 

1. The case d + c < 6 with 2 < d < 4 and 2 < a < 4 

To study its solutions near the bulk critical point, we 

must know how Wd.n behaves for small r. Since the 
long-ranged interaction (x b does not modify the lead- 
ing infrared behavior, it is justified to expand in b. A 
straightforward calculation (see Appendix A) shows that 
provided 2<d<4, 2<cr<4, and d + a <6, 

Wd,o(r|^)-l^d,o(0|cx)) 

-Ad r""!^-^ ^{wd + b Wd.a) r + 0{r^) 



(d+o-)/2-2 



[l+o{r)]+0{b^), (4.11) 



where 



and 



Bd,(T = 



7r((i + cr-2) 
2 (47r)'i/2r(d/2) sin[7r(d + a) /2] 



> . (4.13) 



We insert the above result into the bulk equation of 
state (4.10), keeping only the explicitly shown contribu- 
tions. The resulting equation for the scaled inverse sus- 
ceptibility Tood^^ is expected to take a scaling form. To 



the linear order of our analysis in b and the irrelevant 
scaling fields g„ and g^^, this is the case provided a term 
linear in b is included in g^. Such a contribution is an- 
ticipated on general grounds because in a t/)^ theory with 
coupling constant u, the RG flow of the running variable 
u{i) should be affected by terms linear in 5; technically, 
this may be attributed to the fact that single insertions 
of the long-ranged operator (2.9) require contributions 
linear in b of the countcrterm. 

On the other hand, the scaling field g^ should have 
no contribution of zeroth order in b because, given an 
initial Hamiltonian without long-range interactions (b = 
0), no long-range interaction can be generated under a 
RG transformation. 

In conformity with these ideas, the choices 



9u;{b)=Wa + bWd,a , 9a{b) = b, 



(4.14) 



(up to nonlinear contributions and a redefinition of the 
scales of these fields) turn out to be appropriate. They 
entail that the resulting bulk equation of state scales, so 
that solutions Too to it can be written as 



9t \9h9t ^9u9t i9a9t 



(4.15) 



where the critical exponents 7 = (2 — 77), i^, 77, A = 
{i'/2){d + 2 — f]), ojcr, and u! take the spherical- model 
values 



7SM = 2l^SM 



d-2 ' 



VSM = , Asm = 



d + 2 
2{d-2) 



Wo-,SM = CT - 2 , (4.16) 

and (2.18), respectively. The function TZ^°°^ is given by 

= nt^ (xh) + n\°-^ {xh) + X, 7^S°°) {xu) 

■^o{x^,Xa), (4.17) 

where TZ^f^\xh) is the solution to the asymptotic scaled 
bulk equation of state 



l + xl[Rt\xH)\ " = Aa[n^^\xH 



(00), 



1 (d-2)/2 



, (4.18) 



(4.12) while the remaining two functions are given by 



{d-2)Ad[n'^^\xn)]'"^^\Axl 



(4.19) 



and 



Xh) 



. {d+a+2)/2 



{d-2) Ad [nf\xh)] 



d/2+1 



(4.20) 



■^xl 
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For zero magnetic field, the above findings simpliiy 
considerably, giving 



' oo|/(=0 



1 



9t 



{d- 2) A A Ad 



2ga{b)Bd,cr /gtV 
{d-2)Ad \aJ 



(4.21) 



where again the spherical-model values (4.16) and (2.18) 
must be substituted for the critical exponents 7, u, uj„, 
and uj. 



Logarithmic anomalies and the case d + a = 6 with 
2 < d < 4 



The above results get modified by the appearance of 
logarithmic anomalies when d = Aord + a = &. Our 
ultimate interest is to understand the consequences this 
has for finite-size scaling and the Casimir force in the 
latter case. Since logarithmic anomalies occur already in 
the bulk theory, it will be helpful to clarify their origin 
first in this simpler context. 

That the finite-size behavior gets modified by the pres- 
ence of logarithmic anomalies when d + a = Q was recog- 
nized already in a paper by Chamati and one of us [37]. 
However, no explanation of their cause within the gen- 
eral context of RG theory was given there. Here we wish 
to fill this gap. As wc shall see, despite some similarities 
with the situation at the upper critical dimension d = 4, 
the mechanisms by which they are produced in the case 
d = 4 and d + a = Q with 2 < d < A are different. 

Let us begin by recalling the well understood case d = 
4 [78]. The coefficients Ad and Wd both become singular 
as d 4 (sec, e.g., Rcf. [79] and Appendix A). Although 
Wd is nonuniversal, its pole part at d = 4 (a single pole) 
is universal and equal to that of Ad, so that the sum 
of these two terms in Eq. (4.11) produces a finite rlnr 
contribution in the limit d ^ A. As a consequence, the 
leading thermal singularity of Too takes the form 



1/1=0 



5t/|ln5t| 



(4.22) 



In the framework of RG theory the appearance of log- 
arithmic anomalies means that the Hamiltonian Ti trans- 
forms under a change of momentum scale 11 ^ fi£ into a 
transformed one 'H{{gj{i)},i) whose ^-dependence can- 
not fully be absorbed through scale dependent scaling 
fields (jj (f.) but has an additional explicit dependence on 
£. We follow here the notational conventions of Wegner 
[78, 80, 81]: The gj{i) are nonlinear scaling fields with 
initial values gj{l) = gj and eigenexponents yj; i.e.. 



g^{e) = £-y^g^ 



(4.23) 



We denote their linear counterparts as and let flo 

with 7/0 = d he the special field associated with the vol- 
ume. If the linearized RG operator is diagonal in the 



variables /ij, then these fields usually satisfy fiow equa- 
tions, which to quadratic order can be written as 



de 



where a,j/j = aikj and = [80] . Provided the condi- 
tions 



Vi Vj + Vk 



(4.25) 



are fulfilled, one arrives at an expansion of the form 

ft = 5i + ^ X] 9j9k + ■■■ (4.26) 



with 



^ijk — 



(^ijk 



Dj +yk-yi 



(4.27) 



Similar conditions involving sums of more than two 
eigenexponents, e.g., jji 7^ yj +yk + yh niust hold in order 
that the contributions of third and higher orders have 
a corresponding form with scale-independent expansion 
coefficients. 

When conditions such as Eq. (4.25) are violated so that 
yi equals a sum of other RG eigenvalues, the coefficients 
of the expansion of the linear fields p,i in the nonlinear 
ones gi become scale-dependent, involving logarithms of 
£ or even powers of such logarithms. For example, when 
yi = yj + yk for a single triple k) with aijk ^ 0, then 
bijk gets replaced by [78, 80] 



bijk{() = -aijkia^ 



(4.28) 
i{gt) such that 



Upon making the usual choice £t 
\9t{^t)\ = Ij logarithms of t result. 

Let us first consider the case d + a < 6 with d,a G 
(2,4), and ignore the contributions from all irrelevant 
fields. Then the inequalities (4.25) as well as the condi- 
tion that the linearized RG operator be diagonal at the 
critical fixed point in the (relevant) fields are satisfied. 

The logarithmic singularities one encounters at the up- 
per critical dimension d = A have two sources [78]: (i) 
The exponent /io = d is equal to twice the thermal RG 
eigenexponent yt = 1/j^; (ii) a marginal operator (cj)'^) 
must be taken into account, so that an infinite number 
of eigenexponent inequalities (4.25) and its analogs in- 
volving more than three RG eigenexponents are violated. 
The known consequences arc that the leading thermal 
singularities have logarithmic anomalies which for gen- 
eral values of n consist of nontrivial powers of In gt. 

Next, we turn to the case d -I- cr = 6 with d,aE (2, 4). 
A similarity with the case d = 4 is that the coefiicient 
'Wd,a (which is again universal) has a single pole at d -|- 
(7 = 6 that cancels with the pole of Bd.a such that a 
contribution (x r In r is produced in the limit cr ^ 6 — d of 
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Eq. (4.11). Thus the analog of this equation for d+a = 6 

becomes 



« -Ad'r''^'^-^ + {wd + Wdb)r 

r— >0 

+ bKdrliir + 0{b^) 



with 



Wd 



Kd 



-I- in''^ 



where Kd denotes the conventional factor 

2 



Kd 



5{\q\-l) 



(47r)<i/2r(d/2) 
while w^^Q_^ means the regular part 

2Kd 



> 



lim 

(7 — *6 — ( 



cr + d-6 



+ Wd. 



(4.29) 



(4.30) 



(4.31) 



(4.32) 



of Wd,a at CT = 6 — d. 

Upon substituting the above result into the bulk equa- 
tion of state (4.10), we see that instead of Eqs. (4.15)- 
(4.20) we now have 

Too « g2[n^^\xh)+grT^'^^\xh)x 

X [wd + bwd + bKdln [g]Tl'^^\xh)] }] , (4.33) 
which simplifies to 



'oo|/i=o 



\AdJ \ ^id-2)Ad\AdJ 

Wd + bwd + b - — - In — 
a — z Ad- 



(4.34) 



when h = 0. 

The origin of the logarithmic corrections oc 6 is due 

to the previously mentioned mixing of the 6-indcpcndent 
linear part of the irrelevant scaling field g^ (which we de- 
note as Ha) with ij,a oc b, which led us to conclude that the 
scaling fields (/^ and g^ can be chosen as in Eq. (4.14) up 
to nonlinear contributions. Recalling that /io- is expected 
to contribute to the change of /x^; under RG transforma- 
tions, but /Xct cannot be generated when 2 < cr < 4 if 
the initial Hamiltonian does not involve any long-range 
interactions, one concludes that this translates into flow 
equations of the form 



~ dl ^ V^t^^ + <^w(7 /U(7 + • • 



(4.35) 



with a^cr 7^ 0. As long as d -|- cr < 6, the eigenexponents 
ya and differ. In that case these flow equations yield 



gM) 



a. 



ya- - Vuj 



(4.36) 
(4.37) 



which in turn implies that g^ involves a linear combina- 
tion of fi^ and /icr, in conformity with Eq. (4.14). 

For d + a = 6, the spherical model yields ya = yw = 
d — 4 [cf. Eq. (2.11)]. Owing to this degeneracy, the ex- 
pansion (4.36) gets replaced by 



Mc^(^) = 5t^(^) - aujaga{tj ln£ 



(4.38) 



which in turns leads to the logarithmic temperature 
anomaly in Eq. (4.34). 

The general mechanism we have identified here as pro- 
ducing the logarithmic anomalies in the case d-\-a = 6 is, 
of course, not new; a brief discussion of it may be found 
in Sec. V.E.I of Ref. [80]. 



C. Finite-size scaling form of equation of state 

We now proceed with our analysis of the finite-size 

behavior. To this end we must work out the large-L de- 
pendence of the functions /^Wd,a. and AUd,(i- Expanding 
again to linear order in b gives 

AWd,n{r\L) = ^W^%{r\L) + 6 Aiy« (rji) + O(b') , 

(4.39) 

where the superscripts (0) and (1) on the right-hand side 
indicate respectively the function AWd,Q{r\L) and its 
first derivative with respect to b, taken at 6 = 0. Prom 
Eqs. (3.7) and (4.4) we obtain 



(d) 



OO 



2 cos{qikL) 
r + n{q) 



and 



A<)(rlL)=£ / 

r._i Jq 



(4.40) 



(4.41) 



The q-integrations (cosine transforms) appearing in 
these equations are well-defined as long as r > and 
1/ > 0. In order to obtain the asymptotic behavior of the 
functions (4.40) and (4.41) for r — > 0, we extend the q- 
integrations to the full q-space R'' and make the replace- 
ment fi(q) q^ in their denominators. This amounts to 
the omission of contributions that are regular in r or less 
singular than those retained. The resulting expression 
for the right-hand side of Eq. (4.40) is easily evaluated 
by noting that it is nothing else than the difference be- 
tween the free propagator G^^^'^^ of Eq. (2.30) and its 
bulk counterpart Goo, given by 



Goo{d\r;x) 



id) 



M-i)li 



Kd/2-i{xVr) 

d/2-l 



(4.42) 



with X = \x\. 
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One thus arrives at 



AWgl{r\L)^^J2Gooid\rL^;k) 
Kd-i 



fc=i 

/-2 dp 



L'^-'^ Jo ^yrL^ + p2 eV^^I^ _ i 



(4.43) 



where the second line follows from the first one with the 
aid of the representation 



Goo{d\r;x) 



2(r + gpV2 



(4.44) 



upon interchanging the summation over k with the inte- 
gration over the d — 1 dimensional wave-vector qy con- 
jugate to £C|| = {x2, • . • , Xd), the component of x perpen- 
dicular to ei. 

In order to compute the analogous approximation for 
AWjj^Q{r\L), we proceed as follows. Applying the iden- 
tity 



(r + (jf2)2 



dr 



r q' 



cr-2 



r + 



(4.45) 



to the integrand of Eq. (4.41), we see that the right-hand 
side of this equation is the derivative dr of an expression 
that differs from the right-hand side of Eq. (4.40) merely 
through an extra power of in the integrand. This 
tells us that the roles of the film propagator GL{d\r;x) 
and its bulk counterpart Gcc{d\r-,x) in Eq. (4.43) now 
are taken over by the modified film propagator 



(4.46) 



and its L = oo analog, respectively, which obviously re- 
duce to the former two when a = 2. 
Let us define the function 



QdAy) = |[GL=i(d,(T|j/;0)-Geo(rf,a|y;0)] 



q^-^ cos(gifc) 

y + q'^ 



(4.47) 



Here the second representation follows again by Poisson's 
summation formula (A4). 

In terms of this function, the analog of Eq. (4.43) be- 
comes 



(4.48) 



where the prime indicates a derivative, i.e., Q'd„{y) = 
dQd,a{y)/dy. Furthermore, our result (4.43) for 



Awj Q(r|Z/) can be written as 

A<)(r|L)«L-('^-)-|^Q,.,(rL2) 



(4.49) 



Explicit results for the propagator Goo{dT(y\r]x) and 
the functions Qd,aiy) are derived in Appendix B. As is 
shown there, Goo{d,a\r;x) can be calculated for general 
values of cr G (2, 4) and expressed in terms of generalized 
hypergeometric functions. From these results the asymp- 
totic behavior of the functions Qd.a{y) for large and small 
values of y can be inferred in a straightforward manner 
(see Appendixes B 2). We managed to express Qd,a{y) for 
general values of (rf, a) in terms of elementary and special 
functions up to a series of the form X)°li(-), but have not 
been able to obtain closed-form analytic results for these 
series in general. However, for a variety of special choices 
(d, a) , we succeeded in deriving explicit analytic expres- 
sions for the functions Qd,cr- In particular, all functions 
Qd,cT required for the analysis of the case d = a = 3 of 
nonretarded van-der-Waals interactions in three dimen- 
sions are determined analytically in Appendix B. 

From the above results the finite-size scaling form of 
the equation of state near the bulk critical point follows 
in a straightforward fashion. Let us choose the scaling 
variables t and h in Eq. (2.19) as 

(4.50) 
(4.51) 



h = ;ji^-i/2i(d+2)/2 



gu and go- in accordance with Eq. (4.14), and introduce 
the scaled inverse susceptibility 

fL = tlL^'" = tlL" , (4.52) 

where again the spherical-model values (4.16) were uti- 
lized for the exponents A/p and j/v. 

Upon subtracting from the equation of state (3.21) 
its bulk analog at the critical point and inserting 
Eqs. (4.2), (4.11), (4.48), and (4.49), we obtain for the 
case 2<d<4, 2<cr<4, and d + a <6: 

i « - (h/fL)^ + Adrp-' - 2fl'Qd.2{rL) - 9..fL 



(4.53) 



The result has the expected scaling form. We can solve 
for vl (at least in principle) to determine it as a function 
TZ of the other scaled variables. Hence we have shown, 
to linear order in g^^, g^, and 6, that the inverse suscep- 
tibility rL can be written as 

rL=L-^n{t,h,gu,,ga) (4.54) 

in the appropriate finite-size scaling regime. By analogy 
with the expansion (4.17) made in the bulk case, we write 

TZ{i, h,g^,ga) 
= Tloit, h) + g^ TZu{t, h) + g^ TZaii, h) + o{g^,g„) . 

(4.55) 

Here TZo{t,h) is the solution to Eq. (4.53) with g^ and 
5ct = set to zero. The other two functions are found to 
be given by 



2[7^o(^,fe)] 
Af{t, h) 



(4.56) 
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and 



n^{t,h) = ' . ' ° (4.57) 



with 



h) 



+ 47eo [Qc^,2(7^o) - 7^0 Qd,2(^o)] , (4.58) 

where TZq stands for 7?.o(t, h). 

In the large- 1/ Hmit the foregoing results must reduce 
to our above ones for the bulk, Eqs. (4.15) and (4.17)- 
(4.21). This implies the limiting behavior 

Tl[i{L)XL),9.{L),9a{L)] 

L^oo 

and corresponding relations between the other TZ- 
functions and their bulk counterparts, namely 

TZaim^HL)] i^"n^r\9h9t'") , a = 0,a;,a. 

L — *oo 

(4.60) 

For d + a = 6 with 2 < d < 4 and 2 < cr < 4, a 
logarithmic anomalies appear again in the equation of 

state and its solution. A simple way to obtain these is 
to take the limits cr ^ 6 - d of Eqs. (4.53)-(4.58). This 
yields 



t 



L 



L2 



(4.61) 



and 



TZoit, h) + L-'^^-'^^n^{t, h)[wd + Wdb 
+ 2bQ'd^e-dMt,h)][1lo{t,h)]-' 
+ bKdln[L-^no{i,h)]}'j . (4.62) 

The logarithmic anomalies manifest themselves 
through the contributions that depend explicitly on 
InZ/^ (rather than merely on scaled variables). 



D. Relation between finite-size and bulk inverse 
susceptibility 

The results of the previous section can be combined 
with those for the bulk equation of state to express the 
inverse scaled finite-size susceptibility in terms of its 
bulk counterpart 



(4.63) 



rather than the scaled temperature field t. The relation- 
ship between and f^o will be needed in the next section 



to determine the excess free energy as a function of the 
inverse bulk susceptibility r^o ■ Since the second-moment 
correlation length ^oo of the spherical model is given by 
Too (up to a normalization factor), this gives us ri, and 
the excess free energy expressed in terms of ^oo- 

1. The case d + a <6 with 2 < d < 4 and 2 < a < 4 

We equate the finite-size equation of state (4.53) with 

its analog for foo , 

t ^ ^h^f-J + Adft-'^/' 

- 9u.f^ + 9aBd,ar^^+''-^^^^ , (4.64) 

and substitute for fz, the ansatz 

vl ~ Ro{foo,h) + gaRa{roo,h) + gojRu;{foo,h) . (4.65) 

This yields for Rq = Ro(^(x))^) the equation 

2Ro-^Q..2(Ro) = Ad{Rt'^/'-f(^-'n 



and for the other functions the solutions 

2R 

Roj{roo,h) = (Ro - foo) 



(4.66) 



(4.67) 



and 



K{roo,h) = — 



-4)/2 



Bd,.{f^^+^^ 
_R(d+.-4)/2^^2Q',,,^Ro)], (4.68) 

where N means the function 

N(foo,/i) = 4/^2 R-2 + (d _ 2)Ad R^,"""''/^ 

+ 4Ro'Qd,2(Ro)-40:i,2(Ro). (4.69) 

2. The case d + cr = 6 with 2 < d < 4 

The analog of Eq. (4.64) is given by Eq. (4.61) with 
fi replaced by foo and the terms involving Qd^2 and 
Q'^ 6-d dropped. Owing to the presence of the logarith- 
mic anomaly oc b, the ansatz (4.65) must be modified so 
as to allow for an explicit L-dependence of R^: 



fi ~ Ro{foo,h) + gaK{foo,'h;L) 
+ {wd + Wdb)L^-'^ Rc.(foo,/i) 

Instead of Eq. (4.68), we now have 

2R r r 

R„{f^,kL) = ^|i^d[Roln(Ro/L2) 



(4.70) 



- foo In 



L2 



2Q'd,6-d(Ro)}, (4.71) 
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where the function N continues to be given by Eq. (4.69). 
Likewise, Eqs. (4.66) and (4.67) for Rq and R,^ remain 
vaHd. 

In the case of primary interest, d = 3, these results 
can be augmented by determining the exphcit solution 
to Eq. (4.66) for h = 0. To this end, we substitute the 
result (B19) derived in Appendix B 1 for the function 
Q3,2- Straightforward algebraic manipulations then lead 
to ' 



Ro(^cx)j 0) = 4 arccsch^ [2 exp ( — \/r^/2)] 

"i ^ v'4 + eV^/2)] . (4.72) 



41n'' 



This and the associated scaling function that follows from 
it via Eq. (4.67) are depicted in Figs. 1 and 2), respec- 
tively. Figure 3 shows the function Rcr{f^,0;L) defined 
in Eq. (4.71). 



o 

8 



3 




FIG. 2: Scaling function Raj(roo,0) for d = 3 one obtains by 
inserting Eq. (4.72) into (4.67). 
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FIG. 1: Scaling function Ro(roo,0) for rf = 3, as given 
by Eq. (4.72). The dotted line represents the asymptote 
Ro,as(roo,0) = foo that this function approaches for large val- 
ues of foo in an exponential manner. 

In conjunction with Eqs. (4.67), (4.68), and (4.70), 
the result (4.72) gives us the asymptotic behavior of rj, 
for ft, = in three dimensions including corrections-to- 
scaling, in an explicit analytic form. 
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FIG. 3: Function Ro-(foo,0;I/) for d = 3 and the indicated 
values of L, as obtained by insertion of Eq. (4.72) into (4.71). 



pendix C, giving 



A[/rf,a(0|i) « L-'^[Ag^(rf) + 5.AGM(rf,^) + o(62)] . 

(5.1) 



Here 



V. FINITE-SIZE BEHAVIOR OF FREE 
ENERGY AND CASIMIR FORCE 

We now turn to the computation of the finite-size free 
energy (3.15), beginning again with the case d + a <&. 



A. The case d + a <& with 2 < d < 4 and 2 < < 4 

In order to use Eq. (3.19), we need the _L-depcndent 
part of Ud,a.{^\L). The calculation is performed in Ap- 



Ag^(rf) = -T^-'^'^T{d/2)C{d) 



(5.2) 



and 



A?,^(rf,a) = - 



2"-"^ C{d + a-2)V[{d + a- 2) /2] 
7r'i/2r(l-CT/2) 



(5.3) 

arc the values of the Casimir amplitudes (2.21) and (2.22) 
for our Gaussian model, where C,{d) is the Riemann zeta 
function. 

Upon exploiting the relation (B5) between the deriva- 
tive of Qd+2,2{'r)/r and Qd,2{r)/r derived in Appendix B, 
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one can readily integrate Eq. (4.49) to obtain 



ir\L) 



(5.4) 

Likewise, Aulj^Q{r\L) follows by integration of Eq. (4.48). 
A simple integration by parts yields 

AU^%ir\L) « L-('^+-2) Q,,,(rL2) . (5.5) 

The above results can now be combined in a straight- 
forward fashion to determine the scaled free-energy den- 
sity /ii"*"^. One gets 



(5.6) 



with 



A'^^{d,a)+QdAr)- 



Bd,, 



^{d+a-2)/2 



d + a 



(5.7) 



where /^°) = f^^\K) denotes the smooth background 
term (3.16). As indicated, the function Y in the first 
equation must be taken at the solution fi of the scaled 
equation of state dY{t, f, h, gu;,ga)/df = 0, Eq. (4.53). 

The bulk free-energy density (per volume) /bk follows 
from this in a straightforward manner. As is shown in 
Appendix B, the functions Qd,2{r) and Qd,ajt2{r) behave 
for large values of r as 



.(<i+l)/4 



r^oo 2(27r)('^-l)/2 



e-^[l + 0{r-'/^)] (5.8) 



and 

QdAr) = -A^,^(rf,<7) 
respectively, where 



+ 0(r-2), (5.9) 



2°T[(d + a)/2] 



(5.10) 



according to Eq. (B33). Though not needed here, the 
value of this coefficient appears in our subsequent anal- 
ysis; it is positive for 2 < a < 4 and vanishes both at 
a = 2 and 4. Since the same applies to A^^(d + 2, a), 
the results (5.8) and (5.9) are in conformity with each 
other. 

Hence neither the term oc Qd+2,2 in Eq. (5.6) nor the 
sum of (5d,(7 and A^^(d, cr) contribute in the thermody- 
namic bulk limit. The remaining terms yield 



bk • 



-9a 



2rcx) 

Bd.a 



^ „d/2 ,9i0_ 2 

d °° ^ 4 °° 



,(d+<T-2)/2 



d + a-2 



(5.11) 



The difference of the right-hand sides of Eqs. (5.6) 
and (5.11) gives us the scaled excess free-energy den- 
sity In the result, the scaling function TZ of 
Eq. (4.54) must be substituted for f, and for Too, we 
have the scaling form (4.15) and the relationship (4.59) 
between the scaling functions TZ and T^oo- Obviously, 
the resulting expression for L'^~^f^l^^ therefore complies 
with the scaling form (2.14). 

To derive and describe what this means in terms of ex- 
plicit results for scaling functions, it is advantageous to 
eliminate the temperature field gt in favor of the inverse 
bulk susceptibility r^c (which in the spherical model is 
related to the bulk correlation length ^00 via Too oc ^^^). 
The advantage originates from the explicit results we 
have been able to get for the dependence of on foo- 
Denoting the corresponding analogs of the scaling func- 
tions X,...,X^ in Eqs. (2.14) and (2.20) hy Y, . . . , Y^, 
we write 



(5.12) 



with 



y{foo,h,g^,ga) = Yo{foo,h) + gc^Yu{foo,h) 

+ gaY^{f^,h) + ... . (5.13) 

The above results in conjunction with those of 
Sec. IV C and IV D yield the scaling functions 



_ 47rQrf+2,2[Ro(rcx>,fe)] 
Ro(r<x),/i) 

_ hP' [foo - Ro(foo,fe)]^ 

2f^ Ro{foo,h) 

and 

Y,{f^,h) = AGM(d,a)+Qd,,[Ro(foo,/i)] 



(5.14) 



(5.15) 



+ Bd,a 



-(d+o--2)/2 p(d+o--2)/2 







d + a-2 

Ro(foo, fe) (d+a-4)/2 



(5.16) 



Note that foo is the full inverse bulk susceptibility, which 
itself has corrections to scaling ^ g^ and g^ according 
to Eq. (4.54). Expanding it in powers of g^, and g^ to 
express /|™^ in terms of f |g^=g^=o would produce contri- 
butions linear in g^ and gcr, in addition to those involving 
Yu and Y^. 
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Behavior at bulk criticality 



At the bulk critical point (bcp) T = T^^ao , h 



0, the 



above results reduce to 



^ex,bcp 



-(d-l) 



aSM / 7 ^\ gajb) 



(5.17) 



where ^5^' ^^"^ ^a'c' ^^'^ spherical- model values 

of the amphtudcs (2.21)-(2.23), are given by [82, 83] 



= i"o(0,0) 

_ pd/2 



47r(3d+2,2(Ro,bcp) 



and 



A™(d) ^ y^(o,o) 



\0,bcp 



Ro,bcp/4 > 



(5.18) 
(5.19) 



A'^^id, a) +Qd A^o,bcp) 



Bd,o 



d + (T 



_ p(d-|-<7-2)/2 
2 "^O.bcp 



(5.20) 



Here Ro,bcp = Ro(0,0) is the d-dependent solution to 
Eq. (4.66) at the bulk critical point. 

Thus at bulk criticality, the scaling fields gc^ and go- 
indeed give leading and next-to-leading corrections to 
the familiar first term involving the Casimir amplitude 
(d) of the spherical model with short-range interac- 
tions. 



X.(roo,0) . = 



(27r 



.l-d 



^(5-d)/2 g-2f V2 



and 



foo-oo (d- 2)2^2 °° 

x[l + 0(f-V2)], (5.23) 
Yo{f^,0)^= -D„{d)f^ +0{f^) , (5.24) 



where Dc{d) is the constant introduced in Eq. (5.10). 
Unlike lo(^(x>)0) and l^(foo,0), which decay exponen- 
tially, the scaling function {fx, , 0) decays in an alge- 
braic manner. 

Thus the contribution due to this latter slowly decay- 
ing term governs the large-i behavior of the excess free 
energy for T > Tc.oo whenever the coupling con- 

stant b of the long-range potential does not vanish. One 
has 



/■sing 
J ex 



DM 



j^-(d+a-l) 



-g,{h)L-^''+'^-'^DM(^) \ (5.25) 



Ad' 

where we substituted Eq. (4.21) for Voo to obtain the sec- 
ond line. This strongly contrasts with the asymptotic 
form that applies in the absence of long-range interac- 
tions: 



f 



■sing 

ex I (,=0 



■<^'-)-*(i:) 



(5.26) 



The case d-\-a = 6 with 2 < d < 4 



2. Behavior for T > Tc,oo and h = 

Next, we consider the case T > Tc^oo and h = 0. As 
L OD, the scaled inverse finite-size and bulk suscep- 
tibilities fi and foo both tend towards -l-oo. Hence, to 
obtain the asymptotic large-L behavior, we must study 
the behavior of the functions Rq, R^j, and in the limit 
f'x oo. Clearly, Ra{frx.,h) foo as foo ^ oo. To de- 
termine the asymptotic large- foo behavior of Ro(roo,0), 
we choose foo so large that the function (5d,2(Ro) in 
Eq. (4.66)) can safely be replaced by the first term of its 
asymptotic expansion (5.8). Solving for Ro(foojO) then 
yields 

Too— ><x) [a — Z) Ad 

(5.21) 



x[l + 0(f-V2)], 



where foo now also is to be taken at /i = 0. 

Using this result together with Eqs. (5.8) and (5.9), one 
can derive the large-foo behavior of the scaling functions 
Yq, and 1^ in a straightforward fashion. One obtains 



Proceeding along similar lines as in the foregoing sub- 
section, one can derive the analogs of Eqs. (5.6), (5.11), 
(5.12), and (5.13). They read 

L''-\fL - Lf(°^) 

1 . A. 



^^Lt- — - — f^ - —Qd+2,2{fL) 



d .d/2 47r 
2 ■ " " 2f L d ''^ rL 

+ \{wd + mb)L''-^fl 



+ 9a- 



K. 



d ^2 , rL 



A'eZAd) + Qd,6-d{fL) + n In 



(5.27) 



2ro 

, Wd + bWd 2 

H r 

I 4 ' OO 



Kd 
4 



9a— In Too , 



(5.28) 



and 



-1 j-sing 



Yo{fx,h) + iwd + Wdb)L'^ '^Y^{fx,h) 
+ gaY^{fx,h;L) + ... , (5.29) 
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where 



A'^^id, 6-d) + Qd,6-d[Mroo,h)] 
+ ^1 [r^oo - 2 Ro{foo,h)] foo In ^ 



+ Ro(roo,/i) In 



(5.30) 



while the functions Yq and 1^ remain given by Eqs. (5.14) 
and (5.15), respectively. Owing to the presence of loga- 
rithmic anomalies, the analogs of the scaling functions Y 
and Ycr in Eq. (5.13) have an additional explicit depen- 
dence on L. It should also be remembered that logarith- 
mic anomalies reside also in the temperature-dependence 
of the 6-dcpcndcnt corrections to scaling of foo • 

The scaling functions Yo{foo,0) and i^(foo,0) for the 
three-dimensional case are plotted in Figs. 4 and 5, re- 
spectively. In Fig. 6, the function Y„{foo,h = 0; L) is 



-0.06 - 



-0.12 



-0.18 




FIG. 4: Scaling function Yo{foo,0) for d = 3 (full line). 
The dashed line represents the asymptote (5.22). The value 
Yo{0, 0) = A^{d = 3) which Yo{foo,0) approaches as foo -> 
is known exactly: According to Ref. [82], it is given by 
A|^(3) = -2C(3)/(57r) = -0.15305 . . .. 

displayed for the case d = cr = 3 and some values of L. 
1. Behavior at bulk criticality 




FIG. 5; Scahng function y^(f oo,0) for d = 3 (full lined). The 
dashed line represents the asymptote (5.23). 
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FIG. 6: Function yCT(f oo, 0; L) for d = cr = 3 and the indicated 

values of L (dashed-dottcd and full Hues). The dashed line 
represents the corresponding asymptote (5.24) with d = a = 
3. 



The leading corrections to scaling now results from 
the 6-dependent contribution involving the logarithmic 
anomaly. 



Let us again sec how these results simplify at the bulk 
critical point. From Eqs. (5.29)-(5.30) one easily deduces 
the asymptotic behavior 



J-sing ^ 



R 



^'^ R2 In 

— f^O.bcp ^n — ^ 



0,bcp 



+ AG^(d,6-d)+Qd,6-d(Ro,bcp) 



SM / ,N Wd + Wdb 



+ 



(5.31) 



2. Behavior for T > Tc,oo and h = 

Turning to the case of T > Tg^oo and h = 0, let us again 
consider the asymptotic behavior of f^]^^ for L ^ oo at 
fixed T > Tc^oo- Upon inserting the largc-foo form (5.21) 
of Rq into the result (5.30) for Y^, one sees that the contri- 
bution in curly brackets decays exponentially and hence 
asymptotically negligible compared to the algebraically 
decaying contribution from the sum of the two terms in 
the first line of this equation. This means that the lim- 
iting form (5.24) carries over to the present case, except 
that we must set cr = 6 — d. Since the expressions (5.14) 
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and (5.15) for the scaling functions Yq and — and hence 
their hmiting forms (5.22) and (5.23) continue to hold, 
the results (5.25) and (5.26) for the leading asymptotic 
behavior of /|™^ when 6 or 6 = 0, respectively, also 
remain valid. 



S. The case d = a = 3 

In Fig. 7 the scaled excess free-energy densities (5.29) 
of the three-dimensional case with nonretarded van-der- 
Waals-type interactions (a = 3) and without those are 
compared for the chosen value L = 50 of the slab thick- 
ness L. For simplicity, we have set the nonuniversal con- 
stants Wd and lid to unity. As can clearly be seen from 
the double- logarithmic plot (b) , the asymptotic behavior 
for large i/^oo when 6 is characterized by the asymp- 
tote (5.24) and differs strongly from its counterpart for 
the short-range case 6 = 0. 

In order to illustrate the effect of the explicit depen- 
dence of the scaled excess free-energy density (5.29) for 
nonvanishing interaction constant 6 on L, we display in 
Fig. 8 linear and double-logarithmic plots of fc™^ foi' 
a variety of values of L, including L = oo. For the sake 
of simplicity, we have set the nonuniversal constants Wd 
and Wd to unity. 



C. The Casimir force 

Using the results of the foregoing subsection for the 
excess free energy, the large-scale behavior of the Casimir 
force (2.3) can be derived in a straightforward fashion. 
Depending on whether (i-|-cr<6ord-|-cr = 6, we have 



smg 



(5.32) 



or 



L-^ [So (f^ , h) + {wd + Wdb) L-^ H„ (f^ , h) 

+ g^b) L-'^ E^f oo,k L) + ...] , (5.33) 

where lu and Wo- take their spherical-model values (2.18) 
and (4.16), respectively. The scaling form (5.32) should 
hold more generally for the n- vector model with 2 < d < 
4 even when d + a = 6, as long as u> and u>cr are not 
degenerate. This applies, in particular, to the case d = 
cr = 3 of nonretarded van-der-Waals interactions, albeit 
with the appropriate (different) values of co and LOa, and 
different scaling functions. 

By taking the derivatives of Eqs. (5.12) and (5.29) with 
respect to L, one can express the above functions Sq,. . . , 
in terms of the functions Yq,. . . , Y„. One finds 



So(roo,/i) 



d — 1 — 2fryadf 



-hdf, 

V 



(5.34) 



'E<r{roo,h) = 



d + LJ -1- 2foodf„ 



d + LVrr-l- 2for,df 



(5.35) 
(5.36) 



and 



d + Lj^ - 1- 2f^df, 



A , 



Ldr 



hdr 



(5.37) 



where again the spherical-model values (2.18) and (4.16) 
must be substituted for w, Ua, and A/v [84]. 

Noting that the limit of Ro as foo and h approach the 
bulk critical point exists, 



Ro(foo,/i) = Ro.bcp + o(foo,/i) 



(5.38) 



one sees that the same applies to the scaling functions 

Yi{foo,h): 

YiifooM = Fi(0,0) -|-o(foo,^) , i = Q,uJ,a. 

Hence the terms in Eqs. (5.34) (5.37) involving the 
derivatives with respect to foo and h yield vanishing con- 
tributions as the bulk critical point is approached. Using 
this in conjunction with Eqs. (5.18)-(5.20) and (5.31), 
one finds that the values of these functions at the bulk 
critical point become 



Ho (0,0) = (d-l)ASM(d), 
H^(0,0) = {d + u;-l)Al%{d), 
H,(0,0) = (rf + a;,-l)ASM(rf,^) ^ 



(5.40) 
(5.41) 

(5.42) 



and 



E,{0,0;L) = 3[AG^(d,6-rf)+Qd,6-d(Ro,bcp)] 



2 + 31n^^) . (5.43) 



To obtain the critical Casimir forces in the cases d+a < 
6 and d + a = Q, we must simply substitute the scaling 
functions Hj in Eqs. (5.32) and (5.33), respectively, by 
their above values at the bulk critical point. 

The asymptotic forms of the Casimir force as L ^ oo 
at fixed temperature T > oo and zero magnetic field 
can be inferred in a straightforward fashion from the cor- 
responding results (5.25) and (5.26) for the excess free- 
energy density. Depending on whether a long-range in- 
teraction (X b is present or absent, one has 



^sing 



keT 



(5.44) 
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FIG. 7: Scaled excess free-energy density (5.29) of the three-dimensional spherical model in a L x oo^ slab with periodic 
boundary conditions for L = 50, plotted versus the fiuitc-sizc scaling variable L/(^oc = \phx (a). The solid line corresponds 
to the case a = 3 of nonrotardcd van-dcr-Waals-type interactions with (7o-(&) = h = 2/^; the dashed line shows results for the 
short-range case 6 = for comparison. In (b) the graphs displayed in (a) are plotted in a double-logarithmic manner. In this 
representation the asymptote (5.24) (dotted line) becomes a straight line with the slope —2. The nonuniversal constants wa 
and Wd both have been set to unity. 
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FIG. 8: Scaled excess free-energy densities (5.29) of the three-dimensional spherical model in a L x oo^ slab with periodic 
boundary conditions and nonretarded van-der-Waals-type interactions (cr = 3). The results for various choices of L including 
L = oo are shown as linear (a) and double- logarithmic plots (b). The asymptotes (dotted lines) correspond to the power-law 
behavior (5.24). The nonuniversal constants Wd and Wd both have been set to unity. 



or the exponential decay 



6=0 



(27rL)('^-i)/2 



(5.45) 



Here again the spherical-model values (4.16) must be sub- 
stituted for 7 and v. 

Figure 9 shows a comparison of the Casimir 

forces (5.33) of a oo^ x L slab of thickness L = 50 with and 
without van-der-Waals-type interactions ((7 = 3), where 
we have again set the nonuniversal constants and 
to unity. The double-logarithmic plot (b) again nicely 
demonstrates the approach to the asymptote ~ and 



the qualitatively different behavior in the short-range 
case. 

In Fig. 10 we illustrate how the scaled Casimir force 

for the case with van-der-Waals-typc interactions (cr = 3) 
varies under changes of the slab thickness L. 



VI. SUMMARY AND CONCLUSIONS 

In this paper we have studied the effects of long-range 

interactions whose pair potential decays at large dis- 
tances as x~'^~'' with 2 < cr < 4. Prominent examples 
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FIG. 9: Scaled Casimir force (5.33) of the three-dimensional spherical model in a L x oo'^ slab with periodic boundary conditions 
for L = 50, plotted versus the finite-size scaling variable i/^oo = \/fao (a). The solid line corresponds to the case tr = 3 of 
nonretarded van-der-Waals-type interactions with 5^(6) = 6 = 2/3; the dashed line represents results for the short-range case 
6 = for comparison. In (b) the graphs displayed in (a) are plotted in a double-logarithmic manner. In this representation the 
asymptote (5.44) ~ (dotted line) becomes a straight line with the slope —2. The nonuniversal constants Wd and Wd both 
have been set to unity. 




FIG. 10: Scaled Casimir force (5.33) of the three-dimensional spherical model in a Lxoo'^ slab with periodic boundary conditions 
and nonretarded van-der-Waals-type interactions (ct = 3). The results for various choices of L including L = oo are shown 
as linear (a) and double-logarithmic plots (b). The asymptotes (dotted lines) correspond to the power- law behavior ~ of 
Eqs. (5.24) and (5.44). The nonuniversal constants Wd and Wd both have been set to unity. 



of such interactions are nonretarded and retarded van- 
der-Waals forces. The latter arc ubiquitous in nature; in 
particular, they are present in fluids. 

Application of the phenomenological theory of finite- 
size scaling revealed that such long-range interactions arc 
of the kind termed "subleading long-range interactions" 
[36] and hence should yield corrections to scaling in the 
critical regime near the bulk critical point. 

For systems belonging to the universality class of the 
n-component model in d dimensions, the associated 

corrcction-to-scaling exponent LUa-, given in Eq. (2.11), 
has a larger value than its counterpart w associated 



with the conventional leading corrections-to-scaling of 

ri > 3 dimensional systems with short-range interactions. 
Hence the corrections-to-scaling governed by oja- are next 
to leading. 

However, irrespective of whether u) is smaller or larger 
than cOa, the subleading long-range interactions yield a 
contribution to the Casimir force that decays in a power- 
law fashion as a function of the film thickness L, both at 
and away from the bulk critical temperature T^^b- Since 
the fluctuation-induced Casimir force one has even in the 
absence of these long-range interactions, for T > Tc^ de- 
cays exponentially on the scale of the correlation length. 
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the contribution due to the long-range interactions be- 
comes dominant for sufficiently large L. 

To corroborate these findings we solved a mean spher- 
ical model with such long-range interactions — and hence 
the limit n — > cx) of the corresponding n- vector model — 
exactly. For general values of ct, d € (2, 4), we confirmed 
the anticipated finite-size scaling behavior, and deter- 
mined the scaling functions to first order in the irrelevant 
scaling fields g^i and g^- 

A crucial, though not unexpected, discovery was that 
the scaling field associated with the conventional lead- 
ing corrections to scaling, g^^, depends on the strength b 
of the long-range interactions. This dependence plays a 
role in the mechanism producing the logarithmic anoma- 
lies by which the finite-size scaling behavior of our model 
turned out to be modified when d+a = 6. In these special 
cases — which include, in particular, the physically im- 
portant one of nonretarded van-der-Waals interactions in 
three dimensions — anomalies of this kind showed up in b- 
dependent (leading) corrections to scaling. We were able 
to clarify their origin (see Sec. IV B 2): They are caused 
by the degeneracy lo = bj^ oi the two correction-to-scaling 
exponents in conjunction with the 6-dependence of guj- 

Thus, three-dimensional systems belonging to the uni- 
versality classes of the scalar cj)'^ model and its 0{n) coun- 
terparts with n < oo should not exhibit such logarithmic 
anomalies because their correction-to-scaling exponents 
w and u)a are not degenerate. 

It would be worthwhile to extend the present work in 
a number of different directions. We have focused our 
attention here on the case of temperatures T > Tc^oo- 
An obvious next step is a detailed investigation of the 
model for temperatures below the bulk critical tempera- 
ture Tc^oo- For the spherical model with periodic bound- 
ary conditions considered here, such an extension, which 
we defer to a future publication, is relatively straightfor- 
ward. 

The scaling forms derived in this paper on the basis of 
phenomenological scaling ideas involved nontrivial criti- 
cal indices, such as 77 and 7/;^ = 2 — 77, and the correction- 
to-scaling exponent iv. Although the exact results for the 
spherical model we were able to present are in conformity 
with the predicted more general finite-size scaling forms, 
they neither permit us to corroborate the appearance of 
a nontrivial value of r] nor to verify the n-dependence of 
Lu for the n-vector model. A desirable complementary 
check of the phenomenological predictions that is capa- 
ble of identifying nontrivial values of rj as well as the 
n-dependence of it and other exponents can be made by 
performing a two-loop RG analysis for small e = 4 — 
We have performed such an analysis; its results will be 
published elsewhere [48]. 

Valuable alternative checks of our phenomenological 
predictions should be possible by means of Monte Carlo 
simulations. Although it is quite a challenge to perform 
accurate Monte Carlo simulations of near-critical systems 
with long-range interactions, suitable algorithms were de- 
veloped and demonstrated to be quite efficient recently 



[85-87]. We therefore believe that accurate tests of our 
predictions via such simulations are feasible. 

An obviously important direction for further research 
is the extension of our work to other than periodic 
boundary conditions, namely, those of a kind giving a 
better representation of typical experimental situations. 
Important examples are slabs with Dirichlet boundary 
conditions on both boundary planes, or more generally, 
Robin boundary conditions. Although some aspects of 
our above findings should carry over to such boundary 
conditions — e.g., the form of the scaled variables encoun- 
tered here and the power-law decrease of the Casimir 
away from the bulk critical temperature — it is clear that 
any quantitative comparison between theoretical predic- 
tions and results of a given experiment requires that ap- 
propriate boundary conditions have been chosen in the 
calculations. We leave such extensions to future work. 



Acknowledgments 

D. D. would like to thank H.W. Diehl's group and 
Fachbereich Physik of the Universitat Duisburg-Essen for 
their hospitality at Campus Essen. 

We gratefully acknowledge the financial support of this 
work by the Deutsche Forschungsgemeinschaft via Grant 
No. Di-378/5, the Bulgarian NSF Academy of Sciences 
(436 BUL 113/118/0-1), and the Bulgarian NSF (Project 
F-1402). 



APPENDIX A: FOURIER TRANSFORM OF THE 
INTERACTION POTENTIAL 

In this appendix we wish to derive the small- 
momentum behavior of the Fourier transform (3.4) of 
the pair interaction J{x) introduced in Eri. (3.3). To 
this end we introduce a lattice constant for each of 
the principal directions of the simple hypercubic lattice 
in d dimension we are concerned with. We assume that 
the lattice has an odd number 2A^„ -|- 1 (with Na € N) 
of layers perpendicular to the Xa axis, so that its linear 
extension along the a:;Q-direction is La = (2A^q, + l)aa. 

With these conventions the lattice Fourier trans- 
form (3.4) of the pair interaction (3.3) becomes 

Ni Na d 

where x^^^ = (jaOa) and we have introduced the dimen- 

sionless momentum components = QaO-a- The mo- 
mentum q takes values in the first Brillouin zone, i.e., 

Qa = 2-KVa/La with Z/„ = -7V„, -7V„ + 1, • • • , A^a- 

The contribution proportional to Ji of Eq. (3.3) gives 
the usual result for nearest-neighbor interactions on a 



22 



hypercubic lattice: 



Ji(qr) = 2Ji^cos(ga) (A2) 

a=0 

= 2jJd-lf + ^^Y.\^"»+Om]]- (A3) 

a=l 

To compute the Fourier transform of the remaining 

part of the interaction (3.3), which wc denote as J2{x), 
it is useful to recall Poisson's summation formula [88, 
p. 31] 



oo ^ oc 



(A4) 



m— — oo 



Applying the generalized functions on both sides to a test 
function f{t) whose support is restricted to [—L/2, L/2] 
gives 

E E / -e^'^-^'/^m. (A5) 

Since we are interested in a system of macroscopic lat- 
eral extent, we take the limits — > oo for all a > 1, 
keeping the associated lattice constants > fixed. 
We thus obtain 



(A6) 

where Va = riai=i ^cx is the volume of the unit cell, and 
the integration is over a slab QJz, = [—L/2, L/2] x R'^"^ 
of thickness L = L\. 

The terms with m 7^ reflect the lattice structure of 
the model and give contributions anisotropic in q space 
(as well as isotropic ones). Owing to the restricted inte- 
gration regime, the m = term also yields q-dependent 
contributions (which, however, are small for large L). 
These anisotropics add to those originating from the 
short-range contribution (A2) and produce, in particular, 
a nonzero value of the coefficient 64.1 of the anisotropic 

terms in Eq. (3.7). 

We have emphasized the importance of long-range van- 
der-Waals type interactions for fluids before. Let us 
therefore consider the case of simple isotropic fluids. For 
such systems it is appropriate to take the continuum limit 
tta — > 0. Then the contributions of the m ^ terms van- 
ish by the Riemann-Lebesgue lemma. We have 



VaJ2iq) 



J. 



(cont) 
2,L 



(9) 



df-xJ2{x)e~ 



iq-x 



aJi, 



(A7) 

In order that ^2(9) have a nontrivial continuum limit, 
the coupling constant J2 must be scaled such that J2/va 
approaches a finite value ^2^°"'^ > 0. 



The Fourier transform has an explicit L-dependence 
due to the restriction of the xi-integration to a finite 
interval. However, the deviation from its bulk analog is 
small, unless L is very small: The integration over the 
parallel coordinates yields a function of that varies 
^ x^'^~^ for large Xi. The error resulting from dxi 
therefore decreases as L~'^, i.e., decays ^ N^"^ when ai > 
0. Let us ignore this L-dependence and determine the 

behavior of its bulk counterpart ^2'^'*^ (?) ^'^^ small q. 

The calculation of the latter is straightforward. The 
required angular integral is 



dfld e'" " = {2TTf'^{qxf-'^'^ [qx) . (AS) 
Performing the remaining radial integration gives 



/(cont).. _ .(cont) ^ d/2 / g \ j^^/2(Pog) , . 

■^2,00 (9) -^2 {^J r[(d + <T)/2] ' ^^^^ 

where if (7/2 is a modified Bessel function. From its known 
asymptotic behavior for small values of q one easily de- 
rives the limiting form 



j(cont) 



^0 + A2 {poqf + Ai {poqf 
-A, {poqr + 0{q^+^q^) .(AlO) 



in which 



Ao 



-2(a-2)A2 = 
7r''/^r(a/2) 
T[{d + a)/2]) ' 



-8(4-cr)(cr-2)^4 



(All) 



while Act is given by the right-hand side of Eq. (3.8). 

The ratios —A2/A0 and Aa-/A2 yield the values (3.11) 
and (3.9) of the coefficients 112 and b, respectively, for the 
case of vanishing nearest-neighbor interaction constant 
Ji. 



APPENDIX B: CALCULATION OF THE 
FUNCTIONS QdAy) 

According to Eq. (4.47) the function Qd,a{y) can be 
represented as 



QdAy) 



E 



(d-i) Ad) 



y + q 



(Bl) 



y^''+^-^^/^J2^oo{d,a\l;k^). (B2) 



k=l 



To obtain the second form (Bl), we have utilized the 
property 

Goo(d,a|r;x) =r('^+'^-4)/2Goo(d,a|l;xV^) , (B3) 
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of the bulk propagator. 

The case cr = 2 is special in that the summation and 
integration over qi in Eq. (Bl) can easily be performed 
to reduce Qd^2 to a single integral, namely 



in terms of the modified Bessel function of the first 
kind and the generalized hypergeometric function 
1^2(0;/?, 7; 2), respectively. Their Taylor expansions 
read 



Qd,2{y) = 



dp 



pu. 



[eVW - 1] v^y + P^ 



(B4) 



Integrals of this kind were also encountered in Krech and 
Dietrich's work [31, 32] on the Casimir effect in systems 
with short-range interactions. 

From Eq. (B4) it is not difficult to derive a useful re- 
lation between Qd+2,2 and Qd.2' 



d_ Qd+2,2 (y) 
dy y 



2{y) 



4:Try 



(B5) 



To do this, one simply must interchange the differentia- 
tion of Qd+2,2{y)/y with respect to y with the integra- 
tion over p, replace the y-derivative of the integrand's 
y-dependent part by a derivative with respect to p^, and 
then integrate by parts. 

Returning to the case of general a, we note that the 
representation (B2) has the advantage of linking the 
asymptotic behavior of Qd,(T{y) for large values of y to 
that of Goo {d, cr| 1; a;). In addition, there are some special 
values of (d, a) for which it allows one to derive closed- 
form analytical expressions for Qd,a in a straightforward 
fashion. We therefore begin by computing the bulk prop- 
agator. 



1. Calculation of the propagator Goo(d, cr|r; a;) 

We start from Eq. (4.46) and perform the angular in- 
tegrations using our previous result (A8). This gives 



Gooid,a\l;x) = ., / dq 



(27r)''/2 J, 



1 + q 



2 Jd-2 (qx) . 



(B6) 

The required integral can be evaluated with the aid of 
Eq. (6.565.8) of Ref. [89] or Mathematica [90]. One 
obtains 



Goo{d,a\l;x) 

TT csc[7r((i + (t)/2] 



p(rcg)/ d X \ 



(47r)d/2 

+ (|)"-".^r"(i.2-|,3- 



d + a x'^' 



(B7) 



where qFi^^^ and o^f are regularized generalized hy- 
pergeometric function which can be expressed as 



oF(^''^)(;d/2;xV4) = {x/2)'-'/' I^{x) 



and 



iF2ia;P,j;z) 

mm 



(B8) 
(B9) 



oF(-«'(;d/2;z) = ^ 



3=0 



j\r{j + d/2) 



(BIO) 



and 



iF^r^\l;2- a/2,3 -{d + a)/2;z] 



E 

j=0 



r(i + 2-(7/2)r[i + 3-(rf + a)/2] 



• (Bll) 



For cr = 2, we recover the familiar result for the free 
propagator of systems with short-range interactions: 

Goo(d,2|l;;c) = (27r)-'*/2 x-^'^-^)/^ JsT^ (x) . (B12) 

The latter is known to decay exponentially; the familiar 
asymptotic expansion of the Bessel functions K^, implies 
that 



Goo{d,2\l;x) 

^~(d~l)/2p-x r-m-1 



x^oo 2(27r)(''-i)/2 



j=0 



i!r[(d-l-2j)/2] 



+ 0{x-"') 



(B13) 



When 2 < (J < 4, the Fourier transform of the propaga- 
tor Goo {d, cr| 1; a;) is not regular in q at q = 0. This entails 
that the propagator decays only as an inverse power of x. 
An easy way to obtain its asymptotic expansion for this 
case is to start from Eq. (B6), do a rescaling g — > Q = qx, 
expand the factor {1+Q .t^^)^^ of the resulting integrand 
in powers of x"^, and integrate the series termwise. This 
leads to the asymptotic expansion 



Goo {d,a\l;x) 

2<T-2 



m—1 

E 



T\j + (d + a -2)/2] i-4y 



^^00 T^d/2^d+a-2 Y{l-j-a/2) x^i 



+ 0(a;-2") 



(B14) 



Let us see how the above results can be employed 
to compute the required Qd,a- To treat the three- 
dimensional case, we need Qd,2{y) for d = 3 and d = 5, as 
well as (53,3, and their derivatives. Since Q3 2(2/) involves 
Qi.2{y), we also determine the latter. For these choices 
of d, the result (B12) reduces to 



Goo(l,2|l;ar) = -e- 



Goo(3,2|l;x) = — e- 



(B15) 



(B16) 
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and 



Goo(5,2|l;x) = l±^e- 



(B17) 



respectively. Upon substituting these expressions into 
Eq. (B2), the series can be summed, giving 



2 exp(Vy) - 1 
Q3,2(y)=-£ln(l-e-^) , 

and 

Q5,2(y) = ^[Li3(e-v^) +^/yLi2(e-^/^)] 
where Lip is the polylogarithmic function, 



(B18) 



(B19) 



(B20) 



(B21) 



As can easily be checked, these results (B18)-(B20) are 
in conformity with Eq. (B5). Plots of the functions are 
displayed in Fig. 11. 
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FIG. 11: The functions Q^Mv) (full line). QsAv) (dashed), 
and Qi,2(2/)/10 (dash-dotted), respectively. 

For other choices of a and d — including a = d = 3 — 
the series (B2) cannot in general be summed analytically. 

This suggests that one has to resort to numerical moans. 
To this end a different representation of Qd,a, which we 
are now going to derive from Eq. (Bl), proved to be 
more effective. Remarkably, this representation enabled 
us to derive even a closed-form analytical expression for 

QsMv)- 

Note, first, that the subtracted q-intcgral in Eq. (Bl) 
is the bulk propagator at a; = 0. Both the limit a; ^ of 



Eq. (B7) as well as the explicit calculation of the integral 

Goo(rf, a\y; 0) = | csc[7r {d + a)/2] . 

(B22) 

When d + a > 4, this result involves analytic continua- 
tion in d since the integral is ultraviolet (uv) divergent in 
this case. The same L-independent uv divergences must 
occur in the first term of Eq. (Bl), so that they cancel 
in the difference. We find it most convenient to handle 
uv divergences of this kind, which occur at intermediate 
steps, by means of dimensional regularization. Readers 
preferring to work with a large-momentum cutoff A are 
encouraged to utilize a smooth variant of it, since a sharp 
cutoff is known to give unphysical results in treatments 
of finite-size effects based on the small-momentum form 
of the inverse free propagator, i.e., of fl(q) [36]. 

Next, consider a term of the series J2qi Eq. (Bl). It 
is given by the integral 



/■OO 

IdAli^y) = ^d-i dqii g] 



4-2 



? 4^np('^-2)/2 
2 , „2 



y + <ii + q\\ 



Its calculation for gi = is straightforward, giving 
7d,<.(0,y) = i^d_ij/('^+'^-^)/2 



(B23) 



2cos[7r {d + a)/2] 



(B24) 



It can also be computed in closed form for qi ^ 0; 
the result involves a hypergeometric function 2F1 and 
algebraic functions of y and qi. Rather than working 
with this expression directly, it is more convenient to split 
off appropriate terms containing the uv singularities they 
contribute to the series in Eq. (Bl). Whether and 
what kind of subtractions are necessary depends on the 
values of d and a for which Qd.o- is needed. The Qd,a{y) 
with the largest values of d+a encountered in our analysis 
of the three-dimensional case with a = 3 are Q5^2 and 
Qa^s- Thus the largest value of -|- ct for which Q^^^ 
is required is 7. Using power counting we sec that the 
strongest possible uv singularity of the bulk integral J^"^^ 
in Eq. (Bl) is ^ \d+a--4:^ Hence it is sufficient to subtract 
from the integral (B23) its Taylor expansion to first order 
in y. This ensures that the difference, summed over qi, 
produces a uv finite result. All poles must originate from 
the subtracted terms and cancel with those of the bulk 
contribution in Eq. (Bl). 

Accordingly, we decompose Id,a{Qi,y) as 



(«i , y) = 1^ (91 , 0) ^ + Zd,a (quy), (B25) 
fe=o 

where iji^^^ denotes the kth. derivative of the function 
Id. a with respect to its second argument. Computing the 
integrals of the Taylor coefficients and the remainder Zd,a 



y 
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yields 



5 + 2k-d-a 



) (B26) 



and 



■^d,(7(9i,y) 



2cos[7r(d + CT)/2] 



y 



(<T-2)/2 



(y + ^2)(3-d)/2 



^2|^^|rf+.-7 r[(rf-l)/2] 

y + r(3 - a/2) 



(rcg) 



d— 1 d + a — 5 qf 



(B27) 



where B (a, 6) and 2^1""^^ are the Euler beta function and 
the regularized hypergeometric function 



2F^'''^\a, b; c; z) = 2^1(0, 6; c; z)/r(c) , 



(B28) 



respectively. 

We now substitute the above results into the represen- 
tation (Bl) of Qd.a, utilizing the fact that series YI^^^q 
of pure powers of qi give zeta functions: 



and interchanges the integration over t with the summa- 
tion over j. In this manner the series ^3^3 can be 
computed, and one obtains 



Q3,3(y) 



y y^ ^ 

12 47r2 



In 



+ 



y3/2 

2^ 



- +Im 
4 



Vy 

2tt 



(B32) 



A plot of this function is shown in Fig. 12. 
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qie27rZ 
7^0 



(27r)^ 



(B29) FIG. 12: The function Q:i,3{y)- The in-set is a logaritiunic- 
linear plot of this function, which illustrates the approach 
to the limiting value -A°^(3,3) = -7r^/90 implied by 
Eqs. (5.9) and (5.3). 



Qd,a{y) = y 



^ 00 



+ 



C(5-rf-a) 
72^0 



5 — d— (T 



/d,a(l,0) 



OT^rf-o) (0,1) . 



. (B30) 



Evaluating this expression for (d, a) = (3, 2) and (5, 2) 
with the aid of Mathematica [90] , we have checked that 
the previous results (B18)-(B20) for Qi,2(?y), Q3.2J and 
05,2(2/) are recovered. It can also be utilized to determine 
Q3,s{y) analytically. To this end one rewrites the series 
coefBcient as 



■^3,3(27rj,y) = 



y 



s/y r 

arccos - 



27ri 



47r2j 27r ^^/y + A'K^f 



= Vy f dt 

Jo 



8772 j {t + 47r2j2) 



(B31) 



2. Asymptotic behavior of Qd.aiv) for small and 
large values of y 

The asymptotic behavior of the function Qd,2{y) for 
large values of y readily follows from the representa- 
tion (B2) in conjunction with the asymptotic expansion 
(B13) of the bulk propagator (B12). The result one finds 
for general values of d, 



y 



Qd,2{y) = ^ 

v^oo 2 (27r) 



(<i+l)/4 



-l)/2 



1 + (y-'/^) 



(B33) 

can be verified to be in accordance with the large-j/ be- 
havior of the explicit expressions (B18)-(B20) of these 
functions for d = 1, 3, and 5. 

To determine the large-y behavior of Qd,a{y) with 2 < 
(7 < 4, we insert the asymptotic expansion (B14) into 
Eq. (B2). The summations over k can be performed for 
the expansion coefficients, giving ^-functions. In this way 
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one arrives at the asymptotic expansion 



QdAv) 



1(7-2 



d/2 



m — 1 r 

E 



T[j + {d + a-2)/2] 



X C(rf + CT + 2j - 2) 



r(l-i-a/2) 

(-4)i 



(B34) 



Again, one can employ our explicit result (B32) for Qa^s 
to verify this asymptotic series. Note, that the se- 
ries (B34) truncates when a is even, e.g., when a = 2. 
This ensures the consistency with the exponential de- 
cay (B33) one has for a = 2. 

The asymptotic behavior of Qd,2{y) for small y can be 
conveniently obtained from Eq. (B30) for 1 < < 7. One 
finds that 



Qd,2{y) 



(d-l)/2 



/^r[(3-rf)/2] 
(4^0^72 
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2n 

(rf-3)C(5-d) o 



(27r) 
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y' + 0{y'), (B35) 



provided d ^ 1,3,5. The behavior in the latter cases 
follows by expansion about these values of d. The poles 
that the r function yields for individual terms at such 
odd integer values of d cancel, and logarithms of y emerge 
when d = 3 or 5. The expansions one gets in this manner, 
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(B36) 
(B37) 



and 
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agree with those of the analytic expressions (B18), (B19), 
and (B20). 

The small-y behavior of Qd.a with 2 < cr < 4 can be 
determined from Eq. (B30) in a similar fashion. One 
obtains 
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(B39) 



To deal with the case of Qa^s, one can set cr = 3 and 
expand about d = 3. This gives 



y 



12 8 87r2 V An^J 



+ 0{y% 



(B40) 



where Ce = 0.7772156 ... is the Euler-Mascheroni con- 
stant. The result is consistent with what one obtains 
from the analytic expression (B32) for Qs.s. 



APPENDIX C: CALCULATION OF AUd,u{0\L) 



We start from Eqs. (3.13) and (4.1), take the thermo- 
dynamic limit L|| — > oo, and utilize the continuum ap- 
proximation. Upon transforming the discrete sum over 
the momentum component qi by means of Poisson's sum- 
mation formula (A4), we arrive at 



oo ,-(d) 

AC/d,a(0|L) = J2 cosijqiL) InQ^ . (CI) 
3=1 ■'i 



We substitute for Vlq its small-momentum form (3.7), 
expand the logarithm as 



InOq = In [g2 + 0(g4)] _ 2 + ^^^2^ ^ ^^2) 



drop all suppressed terms, and extend the q-integration 
to W^. 

The contribution from In is known from the short- 
range case [28, 31, 32], easily calculated, and given by the 
6-independent term on the right-hand side of Eq. (5.1). 
The 0{b) contribution involves a difference of critical 
bulk and finite-size propagators at a; = for which one 
obtains, using Eqs. (B7), (BIO), and (Bll), 



^ OO 

^ {Gl - G^){d, a + 2|0; 0) = ^ G^d, a + 2|0; jX) 



2-d-a 2^-^T[id + a-2)/2]ad + a-2) 

77-^/2 r(l - ct/2) 



(C3) 



Adding both contributions yields the result displayed 
in Eq. (5.1). 
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